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1.0 SUMMARY 


The program L217 (EOM) is designed to formulate a set of second order linear 
differential equations that describe the motion of an airplane relative to its level 
equilibrium flight condition. The equations are formulated using the Lagrange energy 
approach and under the assumptions that: 

• The airplane is in straight and level flight. 

• All motions are small. 

The option does exist to replace the structural damping with its equivalent viscous 
damping representation. Inertial coupling effects can also be included in the generalized 
inertia matrix to account for any control surface freedoms that may have been added to 
the basic set of structural modes. 

Aerodynamic data may come from either the FLEXSTAB computer program system or 
the Doublet Lattice L216 aerodynamic program. The FLEXSTAB aerodynamic data can 
only be used with the quasi-steady aerodynamic formulation. The Doublet Lattice 
aerodynamic data can be used with either the quasi-steady or the full unsteady 
formulation. 

Output from L217 (EOM) consists of: 

• Equations of motion coefficients for either the quasi-steady or the full unsteady 
formulation. 

• Aerodynamic forces for each aircraft component, such as wing, fuselage, vertical 
tail, etc. 



Significant program restrictions are: 

e Total number of degrees of freedom cannot be greater than 100. 

• Total number of aerodynamic elements must not exceed 400. 

e Maximum number of interpolation points associated with any one aerodynamic 
surface is 100. 

• Maximum number of gust zones (regions for calculating gradual penetration 
effects) is 35. 

• Maximum number of reduced frequencies is 20. 
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2.0 INTRODUCTION 


The Equations of Motion program (L217) was developed for use as either a standalone 
program or as a module of the program system called DYLOFLEX (ref. 1). L217 (EOM) 
was designed to meet the DYLOFLEX contract requirements as defined in reference 2, 
These requirements specify the need for a computer program that can assemble and 
generate the data needed to formulate the equations of motion for a flight vehicle and to 
provide the necessary aerodynamic forces that can be used in the flight loads analysis of 
that vehicle. 

The objective of this volume is to assist individuals wishing to use this program. To 
meet this objective, the following items are discussed: 

• Engineering and mathematical equations used to formulate the problem. 

• Program structure and design. 

• Guidelines for the execution of the program. 

A sample problem is also presented to aid the user in the execution of the program. 


3.0 SYMBOLS AND ABBREVIATIONS 


The following list contains items that appear throughout the text except for section 6.3 
(card input). 

Engineering 

notation Definition 


[A].[B].[C] 

[A] 


Partitions of the augmented generalized inertia matrix. 
Diagonal matrix of element areas. 


[AIC] 


Aerodynamic influence coefficient matrix. 




Slender body projected area. 


bj. 


Reference semichord. 




Lifting pressure coefficients. 




Lifting pressure coefficients due to the response motion of 
the airplane and gust encounter, respectively. 








Slender body segment pressure coefficients due to response 
motion and gust encounter. 

The magnitude of the slender body segment’s pressure 
coefficient (at the segment evaluation point) due to Y and Z 
doublets, respectively. 

Upper surface pressure coefficients on interference body 
elements. 


{C2}. {C3} 


Generalized forcing function coefficients. 


[CP1r],[CP2r] 


Lifting pressure coefficients for generalized coordinate 
displacements and rates. 


[CPg] 


Lifting pressure coefficients for gust zone velocities. 


[CPR] 


Lifting pressure coefficients for generalized coordinate 
responses. 


[CPG] 


Complex lifting pressure coefficients for gust zone velocities. 


[D] 


Doublet Lattice matrix of normalwash factors for thin and 
interference body boxes. 
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d 82 /dx,d 8 y/dx 


Derivative of surface local z and y displacement in the local 
X direction. 


[F] 


Doublet Lattice matrix of normalwash factors for slender 
body on lifting surface boxes. 


f 


Cyclic frequency (cps). 


{Fa} 


Aerodynamic forces. 


[FAC] 


Diagonal factor matrix that doubles the aerodynamic forces 
for slender bodies located on the plane of symmetry. 



Distances from the gust reference point to the gust zone 
control points. 


[FG] 


Matrix of aerodynamic gust force coefficients for each gust 
zone. 


[FR] 

{F^},{F"} 

[F1r],[F2j^] 

[F12] 

{F'12} 

[GC] 

to 

i 

k 

Ls 

[“csl 


Interference element aerodynamic force normal to the 
element surface. 

Matrix of aerodynamic force coefficients for generalized 
coordinate response. 

Aerodynamic forces in the local y and z direction. 

Aerodynamic force coefficients for generalized coordinate 
displacements and rates. 

Complex matrix of aerodynamic force coefficients for 
generalized coordinate responses. 

Complex matrix of aerodynamic force coefficients due to gust. 
Gust correlation matrix. 

Structural damping coefficients. 

Imaginary number = 

Reduced frequency = wbrA^T (radians). 

Slender body segment length. 

Inertia matrix for control surfaces. 
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[Ml] 

Generalized stiffness matrix. 

[M2] 

Generalized damping matrix. 

[M3] 

Generalized inertia matrix. 

[M4] 

Generalized aerodynamic stiffness matrix. 

[M5] 

Generalized aerodynamic damping matrix. 


Surface unit normal. 

NDOF 

Total number of degrees of freedom. 

NDOFcs 

Number of control surface degrees of freedom. 

NDOFgx 

Number of structural degrees of freedom. 

NSB 

Number of slender bodies. 

NTB 

Number of thin bodies. 

An 

Virtual displacement in the direction of the aerodynamic 
force. 

[PC] 

Interference body /slender body element correlation matrix. 

P/.Pu 

Lower and upper surface pressure. 

AP 

Pressure difference: P/ - Pu- 

q 

Generalized coordinate. 

q 

Dynamic pressure = l/2pVrp. 

{Q} 

Generalized aerodynamic forces. 

[QR] 

Matrix of generalized aerodynamic force coefficients for 
generalized coordinate responses. 

{Q'} 

Complex matrix of generalized aerodynamic gust force 
coefficients. 

Aq 

Virtual displacement of a generalized coordinate. 

[R] 

Euler rotational transformation matrix. 

[rs] 

Matrix of slender body segment radii. 
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[rs'l 

[SF] 

T 

U 

Vf 

VG 

'^Gy/L 

{w} 

{wa} 

{-wb} 

AW 

X, Y.Z 
Xc,Yc.Zc 

*y-y/>v 

Xi.Yi,Z, 

XG 

{XGcp} 

XoA.YoA.ZoA 

Xras.Yras.Zras 

XsH.YsH.ZsH 

«G 


Matrix of slender body segment radius slopes = [drg/dx] . 

Gust rotation matrix for a symmetrical or antisymmetrical 
analysis. 

Matrix of symmetry factors. 

Kinetic energy. 

Strain energy. 

True velocity. 

Gust velocity. 

Magnitude of a sinusodial varying gust (vertical or lateral). 

Normalwash, component of the flow velocity normal to the 
surface. 

Augmented normalwash. 

Normalwash induced on lifting and interference surfaces by 
the slender bodies. 

Virtual work. 

Rigid body linear displacements. 

Axis correction values (reference axis system). 

Local axis system coordinates. 

Inertial axis system coordinates. 

Gust reference point deflned in the reference axis system. 
Gust control points defined in the reference axis system. 
Coordinates of the aerodynamic axis system origin. 

Reference axis system coordinates. 

Aerodynamic axis shift values. 

Gust angle of attach = vgA^T- 
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Circumferential coordinate for slender body pressure 
distribution. 
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Element dehedral angle (Doublet Lattice). 


Rigid body rotational displacements. 


Interference body dihedral angle (FLEXSTAB). 

«tb 

Thin body dihedral angle (FLEXSTAB). 


Rotational displacements in the local axis system. 

P 

Density of air. 

<P(t) 

Wagner function. 


Kussner function. 

(0 

Angular frequency (rad/sec). 


Interference element summation matrices for the y and z 
slender body directions. 

Ss/A 

Gust delta function: 

= 1 for symmetric analysis. 

= -1 for antisymmetric analysis. 

^x»®y>®z 

Linear displacements in the local axis system. 

[«w] 

Sign function relating normalwash to flow incidence. 


Gradual penetration forcing function. 


Mode shapes for displacement normal to the surface. 


Mode shape for surface streamwise slopes = [d0 /dx] . 


[1] 

Unity matrix. 


Flow incidence matrix. 

Subscripts 

Definitions 

Aug 

Augmented. 


CP 


Control point. 



cs 


Control surface. 


E 

FP 

G 

I 

IB 

Im 

L 

I 

n 

O 

P 

R 

RAS 

Re 

S 

S/A 

SB 

ST 

TB 

u 

Ua 

V 


Excitation. 

Force point. 

Gust. 

Inertial axis system. 
Interference body. 

Imaginary part. 

Lateral gust. 

Local axis system. 

Local structural axis system. 
Normal to the surface. 

Off the plane of symmetry. 
Pitch. 

Response quantity. 

Reference axis system. 

Real part. 

On the plane of symmetry. 
Symmetric/antisymmetric. 
Slender body. 

Side translation. 

Thin body. 

Upper surface. 
Unaugmented. 

Vertical gust. 


VT 

Y 

y 

z 

Superscripts 

y 

z 

t 

• 

Matrix 

Symbols 

[1 

[f 

[]-‘ 

[') 

{} 

I j 


Vertical translation. 
Yaw. 

Local y direction. 
Local 7. direction. 

Local y direction. 
Local z direction. 
Differentiation, d/dx 
Differentiation, d/dt. 


Rectangular matrix. 
Transpose matrix. 
Inverse matrix. 
Diagonal matrix. 
Column matrix. 

Row matrix. 
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4,0 ENGINEERING AND MATHEMATICAL DESCRIPTION 


4.1 FORMULATION OF THE EQUATIONS OF MOTION 

The equations describing the motion of an airplane relative to its level equilibrium 
flight condition are formulated using the Lagrange energy equation (ref. 3). 


d / dT\ dT dV . , 

dt \dqi) dqi ^ dqj ^ , . . . . n 

where: 


( 1 ) 


T = Kinetic energy of the system 
U = System’s internal strain energy 
Qi = Set of generalized forces 
qi = Set of generalized coordinates 

The airplane’s motions are specified by the array of generalized coordinates qi, i = 1, 2, 
n. The kinetic and strain energies are calculated in terms of the airplane’s mass and 
stiffness properties and the generalized coordinate velocities and displacements. The 
generalized forces consist of forces arising from damping effects and aerodynamic forces 
resulting from airplane motions and atmospheric turbulence. Placing the restrictions on 
the analysis that: 

• The airplane must initially be in straight and level flight 
% All motions are small 


equation (1) then leads to a set of n second order linear differential equations. These 
equations can be either one of two forms, depending upon the type of representation 
used to describe the aerodynamics of the airplane. 


The first form of the equations of motion comes from the use of the quasi-steady 
assumption in the formulation of the aerodynamic forces. Under this assumption, the 
aerodynamic forces are calculated by using the steady state aerodynamic characteristics 
with the instantaneous angle of incidence (ref. 4). Unsteady aerodynamic effects are 
approximated by convoluting the instantaneous incidence angle with Kussner and 
Wagner indicial lift growth functions. The resulting set of linear differential equations 
written in matrix form is: 


([Mi]{q} + [M2]{q} + [Ma]{q}) + ([M4] {<!}*<*>(« + [M5]{q}*4>(t)) 

= {C3} 


( 2 ) 


where: 


[Mi],[M2],[M3] 


[M4],[Ms] 

{C3} 


= Generalized structural stiffness, damping, and inertia 
matrices, respectively 

= Generalized aerodynamic stiffness and damping, respectively 
= Generalized forcing function coefficients 


<I>(t), iff(t) 


Wagner and Kiissner indicial lift growth functions, 
respectively 


All matrices are constant coefficient matrices; that is, they are not functions of 
frequency. 


The second form of the equations of motion results from the use of a full unsteady 
aerodynamic formulation. This approach yields a set of matrix equations of the form: 

([Mi]{q} + [M2]{q} + [Mg]{q}) + ([M4]{q} + [MgJIq}) = {Cg} «g (3) 


The matrix coefficients of equation (3) have the same definitions as those of equation 
(2). However, the aerodynamic coefficients are now of a different form. The aerodynamic 
stiffness and damping matrices, [M4] and [M5], are frequency dependent, and the 
generalized forcing function matrix, {C3}, is frequency dependent and complex. 

The matrix coefficients of equations (2) and (3) will be calculated for a whole airplane. 
L217 (EOM) can accept structural matrices for either half or whole airplane. Matrices 
for half an airplane will be multiplied by two. 

Detailed discussions of the formulation of the matrix coefficients will appear in 
sections 4.3 through 4.8. First, however, it is necessary to define the axis systems and 
sign conventions assumed by the program, 

4.2 AXIS SYSTEMS AND SIGN CONVENTIONS 

In the Equations of Motion (L217) program, three types of axis systems are used in the 
formulation of the problem. All types are right-handed systems. The basic set of axes is 
the reference axis system (RAS). This arbitrarily described axis system is used to 
describe the airplane configuration. It is important that the same RAS be used 
throughout an analysis. 

The second type of axis system is the local axis system (LAS). Each airplane component 
has a LAS associated with it. All the motion for a component is assumed to be defined 
in its respective local axis system. All aerodynamic calculations are performed in the 
LAS. All local axis systems must have an x-axis lying in the freestream direction with 
the positive x direction pointing aft. For lifting surfaces, the local x-y plane must lie in 
the plane of the surface with the positive y-direction pointing in the direction of 


12 



increasing span. For slender bodies (and their associated interference bodies), the z-axis 
orientation will be a function of the type of motion the slender body will be allowed to 
experience. Slender bodies that undergo vertical motion (known as z-bodies) must have 
their local z-axis parallel to the Z reference axis. Slender bodies that experience lateral 
motion (known as y-bodies) must have their local z-axis parallel to the Y reference axis. 
Figure 1 illustrates typical reference and local axis system relationships. 



Figure 1. — Typical Axis System Orientations 


In the discussions to follow, two t 3 rpes of local axis systems will be encountered - the 
local structural axis system and the local aerodynamic axis system. These two systems 
are a result of the differences between the structural idealization and the aerodynamic 
idealization. A more thorough discussion of the differences between the two types of 
local axis systems will be given when local axis shifting is considered (sec. 4.9). 

The last axis system that enters into the formulation is the inertial axis system. The 
L217 (EOM) program assumes that the interial axis system lies in the plane of 
symmetry of the airplane and is oriented such that x is positive pointing toward the 
nose of the airplane and z is positive pointing downward (see fig. 1). All airplane rigid 
body motion is assumed to be referenced to the inertial axis system. 

The rigid body displacement sign convention is shown in figure 2. Note that with this 
convention positive plunge is down; positive side slip is to the right (facing forward); 
positive pitch is nose up; positive roll is right wing tip down; and positive yaw is right 
wing tip moving aft. Double-headed arrows indicate a right-handed rotation about a 
given axis. 



X 





Figure 2. — Rigid Body Displacement Sign Convention 


The sign convention for the local axis system displacements is shown in figure 3, with 
all surface slopes shown in their positive sense. 



Figure 3. — Local Surface Displacement Sign Convention 
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The relationship between the surface slopes and the surface rotational displacements is: 

d8z 


d7 " 


dx 


(4) 


In dealing with the different aerodynamic formulations and the different types of 
aerodynamic bodies, two types of pressure notations will be encountered. The 
FLEXSTAB aerodynamic approach will give both surface pressures on slender and 
interference bodies and surface pressure differences (lifting pressures) on thin bodies; 
whereas, the Doublet Lattice approach (L216) will always give lifting pressures for all 
bodies. In both approaches, pressure is always compressive (see fig. 4) and the lifting 
pressure is defined as AP = Pj -Pu. 





Figure 4. — Surface Pressures 

4.3 STRUCTURAL MATRICES 

The structurally related matrices of equations (2) and (3) consist of the generalized 
structural stiffness matrix, [Mj], the damping matrix, [M2], and the inertia matrix, 
[M3]. Since the generalized stiffness and inertia matrices are generally formed at the 
time of the structural vibration analysis external to the DYLOFLEX system, no 
calculation of these matrices is performed within L217 (EOM). 

Structural damping is generally assumed to be proportional to displacement, but in 
phase with velocity. This representation is accomplished by making the generalized 
stiffness matrix, [Mj] , complex. The imaginary term of [Mi] is expressed as: 

i [gSD][Mi]{q} (5) 


= Matrix composed of the diagonal elements of the generalized stiffness 
matrix, [Mi] 

= Diagonal matrix of damping factors 


where: 

[»*i] 

[ssd] 


IS 


L217 (EOM) does not calculate an [M 2 ] damping matrix when structural damping is 
chosen. This representation is generally associated with harmonic motion and is built 
into the solution program, L221 (TEV156) (ref. 5) 

The option does exist to transform the structural damping into an equivalent viscous 
damping representation. In this instance, an [M 2 ] matrix is calculated by L217 (EOM) 
as: 

[M 2 ] = [gSD][Mi^/2][M^/2] (6) 

4.4 INERTIAL COUPLING 

In the analysis of an airplane configuration, many flight conditions are usually 
examined. Often the only difference between some of these conditions is that in one 
instance control surfaces may be free to move; while, in another instance, the control 
surfaces remain fixed. The capability exists within L215 (INTERP) (ref. 6) to modify the 
basic set of mode shapes to include the added freedoms of control surface motions about 
their hinge lines. This being the case, the generalized inertia and stiffness matrices 
then become incompatible with the modified modes. Therefore, the option exists within 
the EOM program to augment the generalized inertial, stiffness, and damping (if 
calculated within this program) matrices, making them compatible with the modified 
modes. 


The generalized stiffness matrix (and damping matrix) will be increased by adding rows 
and columns of zeros in order to bring the matrix to the proper size. This is shown in 
equation (7). 


[Ml] 


Aug 


NDOFst NDOF^g 
NDOFst [~[Mi]ua I [0]“ 
NDOFcsL [0] i 


(7) 


where: 


[Ml] 

[Ml] 


Aug 

Ua 


Augmented generalized stiffness 

Original (unaugmented) generalized stiffness matrix obtained from the 
vibration analysis 


NDOFgT = Number of structural (or basic) degrees of freedom that are determined in 
the vibration analysis 


NDOFcs “ Number of control surface freedoms that were added in L215 (INTERP) 

The augmented generalized mass matrix contains nonzero rows and columns that 
represent the inertial coupling effects resulting from the additional freedoms. The 
augmented M3 matrix can be written as: 
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[Ms] 

*■ •*-*Aug 

where: 

[“3)ua 

[A],[B],[C] 


NDOPg^ NDOF(^g 


NDOFsx 

[^S^Ua 

1 [B]1 
1 

NDOFcs 

_[A] 

’i [c]_ 


Augmented generalized inertia matrix 
Unaugmented generalized inertia 
Inertial coupling matrices 


( 8 ) 


The inertial coupling matrices are formed by first assembling a diagonal inertia matrix 
for all mass points affected by the added control surface freedoms. This matrix, 
designated as the control surface inertial matrix, [mcs]> will contain all the needed 
mass and rotational inertial properties for each node. The inertial properties input into 
the [m^s] matrix must meet the following requirements: 


• The inertial matrix can include only those inertial properties that relate to the 
nodal displacements allowed. 

• All inertial data must be for the whole airplane. 

• All rotational inertials must be parallel and/or perpendicular to the freestream. 


The mode shapes describing the displacement at each affected node are obtained by 
using the appropriate interpolation (SA) arrays defined in L215 (INTERP) and 
interpolating for any combination of surface vertical translation, the surface streamwise 
rotation, and the surface rotation perpendicular to the streamwise direction. The 
combination of nodal displacements used for a control surface is left to the discretion of 
the analyst. In submatrix form, the total set of mode shapes for the control surface 
nodes is: 


ui = r<^l \u] (9) 

The ordering of [0]Total corresponds to the order of the inertial data in [mcs]* The 
inertia coupling matrices can now be defined as: 


[A] 

= Wjg [mCs][^]gT 


[B] 

= [A]"' 

(10) 

[c] 

= Wjg [™cs][<^»]cs 

4.5 AERODYNAMIC OPTIONS 



In L217 (EOM), the user has the capability of inputing aerodynamic data generated by 
either one of two aerodynamic methods, FLEXSTAB (ref. 7) or Doublet Lattice L216 
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(ref. 8). Together, these two methods offer the analyst four aerodynamic options in 
forming the equations of motion. 

• FLEXSTAB AIC - quasi-steady approach. 

• Doublet Lattice AIC - quasi-steady approach. 

• Doublet Lattice AIC - full unsteady approach. 

• Doublet Lattice Pressures and Generalized Forces - full unsteady approach. 

In the FLEXSTAB AIC approach, the steady state aerodynamic influence coefficient 
(AIC) matrix is used to formulate the quasi-steady aerodynamics. These quasi-steady 
aerodynamic forces are used to form the generalized aerodynamic stiffness and damping 
and the generalized forcing function of equation (2). Data from the FLEXSTAB program 
system consists of the AIC matrix and the geometry data describing the aerodynamic 
idealization. Also, this approach requires the use of the modal interpolation arrays, SA 
arrays, from L215 (INTERP). 

The Doublet Lattice AIC quasi-steady approach is similiar to that of FLEXSTAB. Here, 
the k = 0 (reduced frequency) aerodynamic data is used to make up the [M4], [M5], 
and {C3} matrices of equation (2). Data from this method consists of AlC-type matrices 
and the needed geometry. This approach also requires the SA arrays from 
L215 (INTERP). 

The third option. Doublet Lattice AIC full unsteady approach, is similar to the previous 
Doublet Lattice option. Here, however, the AlC-type matrices are input at each k value 
for which one wishes to define the equations of motion. This aerodynamic data is used 
to form the nonconstant coefficient matrices of equation (3). This approach also requires 
the geometry data and the SA arrays previously mentioned. 

The last option is that of using the Doublet Lattice-generated pressures and generalized 
forces. It is a full unsteady approach that is used in forming the nonconstant coefficient 
matrices of equation (3). This option requires the geometry data, aerodynamic pressures, 
and generalized forces calculated by Doublet Lattice, but it does not require the SA 
arrays from L215 (INTERP). 

The advantage of using the Doublet Lattice AIC option rather than the Doublet Lattice 
pressures and generalized forces approach is that the former can be used with any 
number of different sets of mode shapes. The AlC-type matrices generated by Doublet 
Lattice are functions of Mach number, k value, and the airplane’s external 
configuration, but are modal independent. If all these variables remain constant, and if 
the only difference between conditions is in the mode shapes used to define the 
structure’s deformation, then the changes in aerodynamic forces due to the changes in 
the airplane’s modal displacements can be calculated without having to rerun the 
complete aerodynamic problem. Since the majority of the cost in solving for the 
aerodynamic forces is the generation of the AIC-t3q)e matrices, this can provide a 
substantial cost savings. 
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The AIC option also has the added advantage of being able to account for gradual 
penetration of the gust. The Doublet Lattice pressure and generalized forces option does 
not have this advantage. The Doublet Lattice program L216, as used in the DYLOFLEX 
system, does not readily lend itself to the inclusion of the gust phase lag effects. 
Calculating gradual penetration effects in L216 (DUBFLX) would require extensive 
modification of the modal input to DUBFLX and the generalized aerodynamic forces 
output from DUBFLX. Therefore, the ability to account for gradual penetration of the 
gust in the DYLOFLEX system was limited to the AIC options of EOM. 

Using the quasi-steady AIC approach, options one and two, offers additional savings in 
that the resulting aerodynamic matrices are also frequency independent. However, the 
user may sacrifice a certain amount of accuracy in predicting the aerodynamic forces in 
some instances where high frequency aerodynamics become significant. 

4.6 FLEXSTAB AIC OPTION 

The FLEXSTAB AIC option uses the geometry data and the steady state AIC matrix 
from FLEXSTAB and the modal interpolation arrays from L215 (INTERP) to calculate 
the aerodynamic forces on the airplane. These forces are then used to generate the 
aerodynamic-related matrices for the constant coefficient form of the equations of 
motion given in equation (2). 

The basic equation in the formulation of the aerodynamics is that of the calculation of 
the pressure distribution over the different aerodynamic surfaces. Using the AIC 
matrix, which relates flow incidence to pressure, the pressure distribution over the 
different elements is given by the following equation: 

{Cp} = [AICJW (11) 


where: 


{Cp} 


= Matrix of pressure coefficients 


[AIC] 


== Steady-state aerodynamic influence coefficients 


{^1 


= Matrix of flow incidences 


The internal ordering and the size of the pressure coefficient and normalwash matrices 
is a function of the type of analysis performed, i.e, symmetrical or antisymmetrical, and 
the type of aerodynamic elements used. These factors determine the AIC matrix 
arrangement. 


4.6.1 INCIDENCE AND NORMALWASH CALCULATIONS 


The normalwash on a surface is defined as the component of the flow velocity normal to 
the surface. Dividing the normalwash by the freestream velocity gives the flow 
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incidence for the surface. The relationship for flow incidence in terms of normalwash is 
defined in FLEXSTAB as: 


^ [^w]{wn} 


( 12 ) 


where: 

A sign function that relates incidence flow and normalwash 
— 1 for all normalwash defined in the local z direction 
-1 for all normalwash defined in the local y direction 

{wn} = Normalwash 

The normalwash caused by the airplane’s motion can be expressed in terms of the 
airplane’s generalized coordinates and mode shapes. 

{wr} = -|vt [</>'sn]cp 


= Mode shapes describing the motion normal to the surface at the surface 
control points 

= Mode shapes describing the surface streamwise slope at the surface 
control points 

= Generalized coordinate displacements and velocities 

= Normalwash due to the response of the aircraft 

The FLEXSTAB AIC matrix requires this normalwash be evaluated at the control 
points of the slender and thin bodies. Normalwash due to interference body motions is 
not allowed by definition of an interference body. 

The mode shapes describing surface normal displacements and streamwise slopes are 
calculated using the SA arrays generated in L215 (INTERP). The interpolated mode 
shapes describe linear displacements in the local z direction of the surface and the 
streamwise slope of these displacements. For thin bodies, this direction is normal to the 
thin body surface. Slender bodies, however, may require motion in either one or both 
the vertical direction (z-bodies) and/or the lateral direction (y-bodies). This is a function 
of the slender body location (on or off the plane of symmetry) and the type of analysis 
performed (symmetric or antisymmetric). Since all SA arrays give motions relative to 
the local z-axis direction, motion for y-bodies must be determined from SA arrays whose 
associated local z-axes are rotated in the lateral direction. The ordering of the elements 
within the normalwash matrix (which corresponds to the AIC column arrangement for 
the symmetrical and anti-symmetrical cases) is shown in figure 5. 


where: 




[■A'Sn] 


CP 


R 


[C] = 
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Symmetric^Ana^^sis 


{w^}, 


Norma Iwash in the slender body local z direction 
for slender bodies on the plane of synmietry. NSBg= 
number of slender bodies on the plane of symmetry 


^” 2 >NSB 


s 


' ^”z/y^SB 




^"z/y^NSBo, 


Normalwash in the slender body local z and y 
directions (the z direction is followed by the y 
for. each element) for slender bodies off the plane 
of symmetry, NSBq = number of slender bodies off 
the plane of symmetry. 


{Wz>i 


{ W } 

^ z^NTB^ 


Normalwash in the thin body local z direction for 
thin bodies off the plane of symmetry. NTBq = 
number of thin bodies off the plane of symmetry. 


Anti*symmetric Analysis 




Normalwash in the slender body local y direction 
for slender bodies on the plane of symmetry 


^''y^SB 
^'"z/y>SB| 
^”z^TB 


{ w } 

^ y^NSB„ 


Normalwash in the slender body local z and y 
directions for slender bodies off the plane of 
symroet ry 


=< 


^"z/y^NSB 


i 


{Wz>i 


Normalwash in the thin body local z direction for 
thin bodies on the plane of symmetry 


^"z^NTBe 

{^J7 


Normalwash in the thin body local z direction for 
thin bodies off tPie plane of symmetry 


5 '^z^ NTBq 


Figure 5. — Ordering of Response Normalwash Matrices for the FLEXSTAB A!C Option 
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The normalwash due to a gust is defined as the component of the gust velocity normal 
to the surface. In L217 (EOM), the gust velocity is defined in the inertial axis system 
with a vertical gust defined positive in the minus z inertial direction (up gust) and a 
lateral gust defined positive in the positive y inertial direction (see fig. 6). 



-► Y, 


X| positive into the page 


Figure 6. — Lateral and Vertical Gust Sign Conventions 

In a symmetric analysis where only a vertical gust is applied to the airplane, the 
component of the gust velocity normal to the surface of a thin body is given by: 


wGtb" 


(14) 


where: 

0TB ~ Dihedral angle of the thin body 

For slender bodies, the local z- direction normalwash is equal to the vertical gust 
velocity,'^Dy , and the local y-direction normalwash is zero. In matrix form, the 
normalwash due to the gust can be written as: 

{wc}g = [R^]g {vGy} (15) 


where: 

[•*ls 

{vGy} 


= Diagonal gust rotation matrix for a symmetrical analysis 

=: Vertical gust velocity at each control point defined in the inertial axis 
system 


In an antisymmetrical analysis, only a lateral gust is applied to the airplane, and the 
component of gust velocity normal to a thin body surface is 


wGtb = ~ ^Gl ®TB 


(16) 
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For slender bodies, the z-direction normalwash is now zero while the y-direction 
normalwash is equal to vq^. In matrix form, the normalwash matrix due to a lateral 
gust can be written as: 

where {G}^ = [Rg]^ {vg^} (17) 


where: 

[RG]^ = Diagonal gust rotation matrix for an antisymmetric analysis 

{vQj^} = Lateral gust velocities at each control point, defined in the inertial axis 
system. 

By referring to figure 5 for the required ordering and using equations (14) and (16), the 
gust rotation matrix shown in figure 7 can be formed for either a symmetrical or an 
antisymmetrical analysis. 

Gradual penetration of the gust can be accounted for by calculating a gust phase lag 
function, The phase lag will be applied to the thin body panels and slender body 
segments in gust zone groupings, not to each individual panel or segment. The gust 
zones are defined by the X-reference coordinate of zone leading edges. Up to 35 zones 
may be used. Each gust zone will have a gust control point; that is, the point at which 
the gust phase lag for that zone will be evaluated. The gust control point is defined by 
the user in the reference axis system. The family of thin body panels and slender body 
segments whose control points fall within the same zone will experience the same gust 
phase lag. The point of zero phase lag is denoted as the gust reference point XG, and 
may be defined anywhere. The distances between the gust zone control points and the 
gust reference point are given in matrix form as: 

{f^} = {XGcp} - XG (18) 

For a sinusoidal gust, the gust velocity at each zone control point is written as: 


f 


r • C/ 1 


'^^V/L 

-1(0^ 
e Vt 


icut 


(19) 


where magnitude of the gust in the vertical or lateral direction. 


To relate the gust at each zone control point to the thin body and slender body 
aerodynamic control points, a gust correlation matrix of ones and zeroes can be formed, 
[GC] . Finally, the gust normalwash for either symmetric or antisymmetric analysis can 
be written as: 


(wg)s^^ - [®*=ls/A 



( 20 ) 
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[BGlj 


(RGl^. 


Symmetrical Analysis 


Ti)i 0 

Slender Bodies on 
the plane of symmetry 

[l/0]i 

• • 
• • 

: : 

Slender Bodies off 
the plane of symmetry 

• 0 • 

. [cos e^g ] 1 ^ 

Thin Bodies off 

• • 

• 

• 

• 

the plane of symmetry 

0 ®TB^NTBo 





ITf] 1 


Anti-symmetrical Analysis 

»••••• 0 



Slender Bodies on 
the plane of symmetry 


Slender Bodies off 
the plane of symmetry 


0 


TB^NTB 


Thin Bodies on and off 
the plane of symmetry 


= total number of 
thin bodies. 

= NTBg + NTBq 


[1^0] = a diagonal matrix whose elements consist of alternating 1's and O's with the first, 
element being a 1. 


Figure 7. — Gust Rotation Matrix for FLEXSTAB 
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4.6.2 AERODYNAMIC FORCES 


The expression for the pressure distribution (eq. 11) can now be used with the 
normalwash defined in equations (13) and (20) to determine the pressure distributions 
due to airplane motion and atmospheric turbulence. The ordering of the elements within 
the pressure coefficient matrix reflects the row order of the AIC matrix. This ordering 
which is the same for both lateral and vertical analyses, is illustrated in figure 8. 

The meaning of the elements of the pressure coefficient matrix differs with each body 
type. Therefore, the calculation of the aerodynamic forces for each body type will be 
discussed separately. 

Slender Body Aerodynamic Forces 

A slender body is a body of revolution divided into a number of streamwise segments 
(see fig. 9). The slender body pressure coefficients are defined at the midpoint of each 
slender body segment. The Cp coefficients represent the magnitude of the 
circumferential pressure variation due to the Y-doublet of the segment at the Y-doublet 
evaluation point. Likewise, the Cp^ coefficient is the magnitude of the circumferential 
pressure variation due to the Z-doublet of the segment at the Z-doublet evaluation point 
(see fig. 10). 

The total circumferential distribution for a segment is the sum of the Y and Z 
distributions. 


Cp(x) — Cp cos/3 + Cp sin^ 


( 21 ) 


where: 


/3 = Circumferential variable used to define the pressure distribution of the Z and Y 

doublets over the surface of the segment. 


Forces in the local y and z directions are obtained by integrating the appropriate 
component of equation (21) over the segment area. 


Fy = “ 5 S Lg Ts So 

Fz = - q St Ts 5^^^Cp(x) cosj3 djS dx 
s 

where: 


( 22 ) 


q 


= Dynamic pressure = 1/2 pVx 


2 


Lg = Length of the segment 


rs 


= Radius of the segment 



fCp*) . 




slender body pressure coefficients due to the Y- 
doublets for all slender bodies, ordering is 
slender bodies on the plane of symmetry followed by 
slender bodies off the plane of symmetry. 



NSB 






Slender body pressure coefficients due to the Z- 
doublets for all slender bodies. Ordering in the 
same as Y -doublet pressure coefficients 


^ P ^NSBI 


SB| 
^SB 




Upper surface interference body panel pressures 
coefficients for all interference bodies. Ordered 
bodies on the plane of symmetry followed by bodies 
off the plane of symmetry. 


{Cp } 

{ACp} 


NM 

1 


Thin body panel lifting pressure coefficients for 
thin bodies on the plane of symmetry followed by 
thin bodies off the plane of symmetry. 


(iCp)i 


(SCp) 


NTBj 


NSB = Total number of 
NIB = Total number of 
NTB = Total number of 


slender bodies 
interference bodies 
thin bodies 


Figure 5 . — Ordering Within the FL EXSTAB Pressure Coefficient Matrix 
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Segment midpoint 


^SB 




Figure 10, — Circumferential Pressure Distribution on a Slender 
Body Segment Due to a Y and Z Doublet 


Performing the circumferential integration, the force expressions reduce to: 

Fy = - q )dx 

s 

Q )dx 

s 


( 23 ) 
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Y Z 

Making the assumption that Cp and Cp ^ are constant over the length of the segment, 
and using the radius value at the segment midpoint, the final expressions for the y and 
z forces for any one segment are: 


Fy— q L§ 7T Cp 

2 

~ q r§ L§ 7T Cp 


(24) 


Defining the slender body segment area as: 


As 2 rg Lg 


(25) 


th 

the y and z segment forces for the J " slender body can be written in matrix form as: 

(Fy) - q- [^%] [A^] {Cp’') 

^ (26) 

{Fzlj - q [-*2J [4]j (Cp^)j 

Interference Body Aerodynamic Forces 

The interference panel pressure coefficients represent the upper surface pressure on 
each panel (see fig. 11). These pressures are evaluated at the panel centroid. By 
definition, interference bodies do not generate their own lift, but are present to account 
for the interference flow about a slender body caused by thin and other slender bodies. 
As a result, the lower surface panel pressure is zero. Thus, the lifting pressure, 
coefficient for each panel is: 



Figure 11. — Interference Body Panel Orientation Angie 
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The FLEXSTAB aerodynamic theory assumes constant pressure panels. Under this 
assumption, the panel normal force is given by: 


Fnjg ~ “ 5 Ajb Cpy 


( 28 ) 


where: 


Ajb = Interference panel area 

£ •' 

Using the orientation angle of the local panel normal with respect to the local z-axis 
(see fig. 11), the y and z components for each interference body panel are: 


9 Aib Cpy 

Fzjg = - g cos fliB AiBCpy 


( 29 ) 


For the interference body, equation (29) for each panel can be written in matrix 
form as: 


o o 

^ ^IbIk [Aib]k {^Pu^K 

Addition of Interference and Slender Body Aerodynamic Forces 


( 30 ) 


To obtain the total segment forces on any slender body that has an associated 
interference body, the interference and slender body forces must be added together. This 
addition is made under the assumption that no interference panels overlap slender body 
segment boundaries; that is, the circumferential interference body segments must lie 
within the corresponding slender body segment. This requirement is illustrated in 
figure 12. Note that the last two circumferential segments of the interference body lie 
within the same slender body segment. 

^SB ^IB 



Figure 12, — Slender Body- Interference Body Alignment 



The addition of the interference panel forces to the proper slender body segments can be 
accomplished by premultiplying the force expressions of equation (30) by a panel 
correlation matrix. This matrix is composed of ones and zeroes and is a rectangular 
matrix of size (number of slender body segments) by (number of interference body 
panels). 




(31) 


where: 

[PC]i 


= S^m of the interference body panel y and z forces at the 
slender body segments 

= Panel correlation that adds the interference body 

panels to the proper J slender body segments 


Combining equations (31) and (30) and defining y and z force summation matrices as: 

[SyJn = [-sin'eiB]^ 

, (32) 

[Sz]k = [PC];^[-cos°0ib]k 

the interference body forces summed to the respective slender body segment midpoints 
can then be written as: 


^ 


(33) 


The total body force now becomes the sum of equations (33) and (26). Referring to the 
Cp matrix ordering shown in figure 7, the total slender body force distribution can be 
written in matrix form as: 


{Fy} 

{Fz} 


SB 


SB 


= q 


[M [ 0 ] [ 2 y] 

[0] [-^ 12 ] [2,] 


[Asb] 


[As“b] 


[Aib]J 







,(CP„).B , 


(34) 


Thin Body Forces 


The aerodynamic forces on each thin body ®^® found by multiplying the lifting 

pressure on each panel by its area. For ths F* ’thin body, this is written in matrix form 


as: 


{Fz}j = q [AtbJj {ACp}j 


(35) 
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Thin body forces act in the local z axis of the respective surfaces. 



i I t 


Total Aerodynamic Forces 


Combining equation (34) with equation (35) for all the thin bodies, the total set of 
aerodynamic forces can be written as: 



[[ 


= q 




[ 0 ] [-^/ 2 ] [ 1 ^] [ 0 ] 
[0] [0] [0] [i] 


[a|b] 


[a|b] 


[a?b] 


[Atb]! 


-1 




■j 


Jj 

{ACpI^B 

^ J 


(36) 


Since the aerodynamic modeling in FLEXSTAB is for half an airplane, the slender body 
segment areas for bodies on the plane of symmetry will be half of the total area. Also, 
for interference bodies on the plane of symmetry, only half the interference shell is 
modeled. Therefore, to obtain whole body forces, the Fy and elements for slender 
bodies on the plane of symmetry are multiplied by two. The final expression for the 
aerodynamic forces thus becomes: 

{Fa} = q [FAC] [2] [A] {Cp} (37) 


where: 

[FAC] = Diagonal factor matrix that doubles the aerodynamic forces for slender bodies 
on the plane of symmetry 

4.6.3 GENERALIZED AERODYNAMIC FORCE COEFFICIENTS 

The generalized forces shown in the Laplace energy equation (1) are formed by the use 
of the principle of virtual work. 


Qi = 


aAW 

3Aqi 


i = 1, . . . , n 


(38) 


where: 

Qj = Generalized force in the i*^^ freedom 
AW = Virtual work 

th 

Aqi = Virtual displacement of the i*^” generalized coordinate 

Virtual work is calculated by multiplying the external forces acting on the system by 
the virtual displacement of the structure at the point of application of the force and in 
the direction of the force. 


. ^ 
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(39) 


AW = LAnjpp{F} 

Expressing the virtual displacement of the structure in terms of the structure’s mode 
shapes and virtual displacement of the generalized coordinates, and performing the 
differentiation of equation (38), the expression for the generalized forces becomes: 

{Q} = Kn]FpW (40) 

In the above expression, the modal displacements, describe the displacement of the 
structure at the point of application of the force (force point) and in the direction of the 
force. 

Generalized Response Aerodynamic Forces 

Using the pressure coefficient equation (11) and the expression for the normal wash due 
to airplane motion (13), the response pressure distribution is given by; 

{Cp^ } = - [AIC] [<A'6 „]cp ^ 

(41) 

= [CP1r] {q}+^[CP2ft] {q} 

Substituting this expression into the aerodynamic force equation (37), the response 
aerodynamic forces are: 

{Fr} = q [FAC] [ 1 ] [A] |[CP1r] {q} + ^ [CP2 r] {q}j 

(42) 

= [FIr] {q} + [F2r] {q| 

Finally, the force expression (eq. 42) is substituted into the generalized force 
equation (40) to give: 

{Qr} = 

where: 

Q 

[SF] = Symmetry factor matrix that doubles the forces of the bodies off the plane of 
symmetry in order to calculate total airplane forces 

It is important to note that the displacements at the force points are defined in the local 
aerodynamic axis system of each surface (see sec. 4.2). Thus a positive displacement acts 
in the same direction as a positive aerodynamic force, resulting in positive virtual work. 


32 



Bringing the response generalized forces to the left side of equation ( 1 ) and writing 
them as coefficients of generalized coordinate displacements and velocities, the M 4 and 
Ms matrices become: 


[M 4 ] = - [SF] [FIr] 

" (44) 

[M 5 ] = - [<t>sXp 

Generalized Excitation Aerodynamic Forces 

The generalized excitation aerodynamic forces are calculated in the same manner as the 
response forces. Using the gust normalwash given by equation (20), the gust pressure 
distribution is defined as: 


{Cpj.) - ^ [AIC] [GC] 


r .1(0-^- 


Vt — iwt 
e yOr 6 

V/L 


' [CPg) 


• u 

1 r.\ 

Vt 


-ICO 




io>t 


and the excitation forces become: 




- 10 ) 


{Fo } = ^ [fAc] [X] [A] [CPg] I e e‘“‘ 


- icot 


- [fg] 

Finally, the generalized excitation forces are: 


(45) 


(46) 


{Qe} = [hXp 



io>t 

e 


= [?] "''Tj vGy/j^ e' 

- {C 3 } e‘“‘ 


4.7 DOUBLET LATTICE AIC OPTION 


(47) 


Basically, the Doublet Lattice AIC option performs the same type of aerodynamic 
calculation as the FLEXSTAB AIC option discussed in section 4.6, (relating the pressure 
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distribution over the aircraft to the flow incidence at several aerodynamic control 
points). However, the differences between the two approaches lie in the method used to 
formulate this relationship and the fact that the Doublet Lattice approach can be used 
for unsteady as well as quasi-steady aerodynamics. 

Unlike the FLEXSTAB approach, in which the AIC matrix contains all the interaction 
effects of slender, interference, and thin bodies, the Doublet Lattice formulation is 
separated into three distinct steps: 

• First, slender body pressures are determined, 

• Second, the normal wash induced on thin and interference bodies by the slender 
bodies is calculated. 

• Finally, the thin and interference body pressures are calculated. 

The Doublet Lattice AIC option uses the modal interpolation arrays from L215 
(INTERP) and the geometry and aerodynamic data from L216 (DUBFLX). The 
aerodynamic data consists of frequency dependent matricies, [F] and [D], which relate 
flow incidence and lifting pressure. 

In the discussion to follow, the full unsteady aerodynamic case will be examined first. 
The quasi-steady formulation will be presented afterward since it is a subset of the full 
unsteady approach. 

4.7.1 SLENDER BODY AERODYNAMIC FORCES 


The Doublet Lattice program does not pass slender body pressures to the EOM program 
for this option. These pressures are calculated within the EOM program using Miles' 
slender body theory employed in L216 (DUBFLX) (ref. 8). Based on this theory, the 
lifting pressure over the slender bodies is given by: 


{ACp} 


SB 


where: 





+ ^[h] 



(48) 




= Array of slender body segment lifting pressures in the direction of the 
normalwash, w^ (n is either local y or z direction) 


[rs].[rs] 


= Array of segment radii and streamwise slope of radii evaluated at 
segment midpoints (r'§ = drg/dx) 


{wn}j{wn} = Normalwash and the streamwise derivative of the normalwash 
evaluated at the segment midpoints. 


Slender bodies that undergo vertical motion (z-bodies) give lifting pressures in the 
z-direction. Y-bodies are slender bodies that experience lateral motion and give lifting 
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pressures in the y-direction. As discussed in section 4.6.1, different SA arrays must be 
used with y-bodies and z-bodies. The ordering within the pressure matrix is all lifting 
pressures in the y-direction (Y-bodies) followed by all lifting pressures in the z-direction 
(Z-bodies). 


r ^ 



ACpY] 

ACn 



J 

— — -V 

P 

I 

SB 




Response Aerodynamic Forces 


SB y 


(49) 


The response pressures on the slender bodies can be calculated by using equation (48) 
with the response normalwash defined in equation (13). The normalwash is evaluated at 
the segment midpoints and is formulated in the same order as the pressure stacking 
shown in equation (49). Making this substitution, the slender body lifting pressures are; 


{ACp } 

R SB 


- f KJcpW) 



(50) 


For simple harmonic motion, velocity and displacement are related by the expression: 


q = io)q (51) 

Using this relationship, equation (50) can be written in complex form as: 

- [[CPKsb]^^ + i [CPKsB],„]{q) 

The slender body forces are obtained by multiplying the pressures by the dynamic 
pressure and the segment projected areas. 


W]J 


q [Asb] + i [cPRsB]jjj,]{q} 

[[FRsbIj^^ + i [FRsbIj^ {q}] 


(53) 
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Excitation Aerodynamic Forces 

The slender body lifting pressures due to the gust are calculated using the gust 
normalwash given in equation (20) with a modified version of the slender body pressure 
equation (48). In the modification of the pressure equation, terms dealing with the 
normalwash time rate of change and streamwise slope are neglected. At high 
frequencies, the unsteadiness produced by these terms violates the assumptions of the 
slender body theory. Therefore, the slender body excitation pressure will be 
approximated by: 




f -“f I 




io>t 


( -iJf. 

lie ^tJv, 


= [CPGsb] 

The slender body gust forces become: 


^Gv/l ® 


i(ut 


{FgIsb = q [Asb] [CPGsb] 
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e Vt 


'^Gv/L ® 


io)t 


= [fGsb] 


- 10 >: 


u. 


Jv, 


^Gv/L ® 


i(ut 


(54) 
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4.7.2 INDUCED NORMALWASH 


The second step in the formulation of the aerodynamics is the calculation of the 
normalwash induced by the slender bodies on the thin and interference body boxes. 
Note that the term "box” in the Doublet Lattice option refers to the smallest division of 
a thin or interference body. It is equivalent to the term "panel” used in FLEXSTAB. 

The Doublet Lattice program, L216, calculates a complex matrix of normalwash factors, 
[F] , at each frequency of interest. The matrix is rectangular, with each row relating to 
a thin or interference body box and each column relating to a slender body segment. 
The induced normalwash at the box control points is given by the product of the 
normalwash factors and the slender body pressures. For the induced normalwash due to 
the airplane response motions, the equation is: 

and due to the gust excitation, the expression is: 


(56) 


(57) 
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4.7.3 THIN AND INTERFERENCE BODY AERODYNAMIC FORCES 


The relationship between the flow incidence at the thin and interference body box 
control points and the box lifting pressures is given by the equation: 


where: 

[D] 


{wa}tb’{^a}jb 


1 

Vt 


{wa} 


TB 


{wa} 


IB 


= [D] 


{ACp} 

{ACp} 


TB 


IB J 


(58) 


= A complex matrix of normalwash factors relating box 
lifting pressures to flow incidence at the box control 
points 

= Augmented normalwashes at thin and interference body 
boxes 


{ACp}^^,{ACp}jg = Thin and interference body lifting pressures. 

The matrix, [D], is a frequency-dependent matrix and is calculated at each frequency of 
interest. The augmented normalwash of equation (58) is the sum of the normalwash due 
to the motion of the bodies or due to the gust encounter plus the normalwash induced by 
the slender bodies. 


Response Aerodynamic Forces 

By the definition of an interference body, that is, a body that does not generate its own 
lift, the normalwash due to interference body motion is set to zero. The thin body 
normalwash due to the airplane motions is given by equation (13). Using the induced 
normalwash expression of equation (56), the augmented response normalwash becomes: 
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'{wrI^b' 

Vx 


“ Vx 

(0 } 


+ [F] {ACp } 

R SB 


[0] 


+ [F] [CPRsb] 


'IB 


(59) 


{q} 


37 


Substituting equation (59) into (58) and solving for the box pressures gives: 


{ACp } ^ 

R TB* 




TB Sz rc3' 


[ 0 ] 


IB 


+ [F] [CPRsb] 


{q} 


= [[CPRTB/IB]j^g + i [CPRtB/Ib]^ J{q} 


( 60 ) 


Im-* 


Multiplying equation (60) by the matrix of box areas and the dynamic pressure gives 
the response aerodynamic forces for thin and interference bodies: 


{Cpj^Itb 

{CpjjliB , 

Excitation Aerodynamic Forces 


{^rItb/ib " ^ t^TB/m] 


(61) 


The augmented gust normalwash is obtained by adding the gust normalwash defined for 
the thin body boxes and given by equation (20) to the induced normalwash due to the 
gust slender body forces given by equation (57). 


{"Aq) 1 

I JtB 


^ [RGtb] [GC] 

Vt ■'s/a ■'tb 
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+ [F][CPGsb] 
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(62) 


Note that by the definition of an interference body, the normalwash on the interference 
body boxes due to a gust encounter is set equal to zero. 

Substituting equation (60) into (58) and solving for the lifting pressures, results in: 
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-1 


n 


[RGtbI [GC]r.,„ 
Vt S/A*- ■'TB 


[ 0 ] 


IB 


= [[CPGtb/Ib]_ + ^[CPGtb/Ib] 
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+ [F][CPGsb] 
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i<ut 

( 63 ) 


Im-* 
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The gust forces on the boxes become: 




{fo)TB/iB “ ^ [■'TB/Ib] [CPGtb/Ib] e ’’ 
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= [FGtb/Ib] 
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4.7.4 GENERALIZED AERODYNAMIC FORCE COEFFICIENTS 


( 64 ) 


Total Aerodynamic Forces 

The total aerodynamic response and gust forces are obtained by summing the 
interference body box forces and the respective slender body segment forces, and then 
multiplying the resultant slender body forces for bodies on the plane of symmetry by 
two. The multiplication is needed since the modeling in the Doublet Lattice program is 
for half an airplane. 


The total response aerodynamic forces can be written as: 


{Fr} = [FAC] 


"[ 0 ] - 

[fRsb] 


+ [2] [FRtb/Ib] 


{q} 




(65) 


Again, using the harmonic relationship between velocity and displacement, the response 
force of equation (65) can be written in terms of coefficients of q and q. Thus; 


{Fr) = [FRlBe W 4 W) 


( 66 ) 


Similiarly, the total aerodynamic gust forces are: 
{Fg} = [FAC] 


[ 0 ] 1 

I + [S][FGtb/Ib] 

[[FGsb] 


. f / 


'^Gv/L ® 


ioit 


(67) 


f 


^Gv/L 


icot 
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The summing matrix in equations (65) and (67) is similar to that used in the 
FLEXSTAB AIC option, but with an arrangement corresponding to the Doublet Lattice 
ordering. In submatrix form, it is shown below: 



[0] [0 ] [0] n 


[0] [Sy] [0] 


]o] [0 ] [s,L 

where: 


[^y] = [PG] [single] 

[2J = [PC] [( 

:os0ib] 


The panel correlation matrix, [PC], performs the same function of adding all the 
interference boxes to the proper slender body segments. In terms of the Doublet Lattice 
box dihedral angle, y, the box unit normal orientation is: 



Generalized Forces 


(69) 


Forming a matrix of normal displacements at the segment and box force points, the 
generalized response forces are given by: 


(Qr) - [SF] (FR] {q} 

- [QR] |q) 


(70) 


and the generalized excitation forces are given by: 


Nr) - [*1 [FG] 




="n"TB/SB 

f 

= m e Vt 


^Gv/L® 


ict>t 


""Gv/L ® 


i(ut 


(71) 
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Bringing the generalized response forces to the left side of the equations of motion and 
writing them as coefficients of q and q, the M4 and M5 matrices become: 


[M4] - - [QR]^^ 

[Ms] = -i[QR],„ 

The C3 matrix is defined as: 

£ 

(cs) - m 

4 . 7.5 QUASI-STEADY FORMULATION 


(72) 


(73) 


The quasi'Steady formulation in the Doublet Lattice AlC option is accomplished using 
the aerodynamic data ([D] and [F] matrices) calculated at a> = 0, which makes them 
real matrices. Also, the harmonic motion relationship between q and q is not used, and 
the normalwash is written as real coefficients of generalized coordinate displacements 
and velocities. 


Setting o) = 0 in equation (50) and grouping terms as coefficients of q and 4, the 
response slender body lifting pressures are given by: 

-rji'»’t*«nlcp*liyt*«„lcp](« (74) 

“ICPlRlsBW-LCPZglgBW 

The expression for the slender body lifting pressures due to gust is given by equation 
(54). 


The augmented normalwash on the thin and interference boxes becomes: 

r-[^. 

Ll 

re. 1 


1 

Vt 




A [F] [CP1„] 


{q} 


(75) 


LL [0] 


IB 


B 

+ 


[F] [CP2^] SB 


(« 


where the [F] matrix is now real. 
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Using equation (75) with equation (60), the response lifting pressures on thin and 
interference boxes are: 

{ACpj^} 1 

(iC r ■ ^ ['^“■Jtbab«> 

With the lifting pressures established, the aerodjuiamic forces and generalized forces 
are obtained by substituting equations (74) and (76) into the proper force equation given 
in section 4.7-4. It should be noted that the gust forces computed by the two approaches 
differ in that the [D] and [F] matrices are real for the quasi-steady assumption. 

4.8 DOUBLET LATTICE GENERAUZED FORCE AND PRESSURE OPTION 

The last aerodynamic option is the Doublet Lattice generalized force and pressure 
option. In this option, the lifting pressures on all slender body segments, thin and 
interference body boxes, and the generalized response forces are calculated by the 
Doublet Lattice program, L216. The L217 (EOM) program uses the lifting pressures to 
calculate the aerodynamic forces on each aircraft component, places the generalized 
forces in the format required for the equations of motion, and uses the response forces to 
generate gust excitation forces. No modal data is required from L215 (INTERP) and the 
formulation of the excitation aerodynamics does not permit gradual penetration effects 
to be taken into account. 

Aerodynamic Forces 

The complete lifting pressure matrix read from the Doublet Lattice program is a 
frequency-dependent complex matrix containing all slender body segment, thin and 
interference body box pressures. The ordering follows that used in the Doublet Lattice 
program, which is: 

• Thin bodies on the plane of symmetry followed by thin bodies off the plane of 
symmetry. 

• Interference bodies on the plane of symmetry followed by those off the plane of 
symmetry. 

• All y-slender bodies on the plane of symmetry followed by all those off the plane of 
symmetry. 

• All z-slender bodies on the plane of symmetry followed by those off the plane of 
symmetry. 
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[CPRtb] 


{^Cp^} 


[CPRib] 

[cpr^sb] 

[cpr^sb] 


(77) 


Multiplying the pressure matrix by the box and segment areas, summing the 
interference body boxes onto the respective slender body segments, and multiplying 
them by the factor matrix which accounts for slender bodies on the plane of symmetry, 
the total set of response forces is: 


{Fr} = q [FAC] [2] [A] [CPR] {q} 

= q [F12] {q} 


(78) 


Using the displacement-velocity harmonic relationship shown in equation (51), the 
forces can be written as coefficients of q and q. 

{FIr} = q[Fi2] {q} 

(79) 

{F2r} 

The gust forces are obtained by developing the relationship between the gust velocity 
and the rigid body displacements and then extracting the proper rigid body columns 
from the force matrix to form the aerodynamic forces due to gust. The sign convention 
shown in figure 6 will be used along with that shown in figure 2 to establish the 
relationships. 


For a symmetrical analysis, the airplane experiences the rigid body motions of plunge 
and pitch. The plunge of the airplane will give the same aerodynamic effect as a 
positive gust encounter. Therefore, the relationship between the vertical gust velocity 
and the vertical translation generalized coordinate is given by: 


VGy = qvT "" 

The pitch generalized coordinate is related to the vertical gust velocity by the equation: 

VGy = Vt qp (81) 


For an antisymmetric analysis, the rigid body motions of side translation and yaw are 
used to derive the gust relationships. The side translation generalized coordinate/lateral 
gust relationship is: 


“ ^ST “ ‘"‘1st 


(82) 
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and the yaw/lateral gust equation is: 


= Vt ( 83 ) 

Equations ( 81 ) and ( 83 ) are used in obtaining zero frequency gust forces, while 
equations ( 80 ) and ( 82 ) are used for all nonzero frequency gust forces. 

Extracting the proper rigid body column from the [F12] matrix shown in equation ( 78 ), 
the gust forces at o> = 0 are: 

{FgJ^.a - 5 {’"' 2 I Vy ■ ''°v/L 


and for > 0 the gust forces are: 

^ \t/ST ^ ^{F'12} VGv/l 


where: 

{F'12} = A complex column matrix of aerodynamic force coefficients formed by 

extracting the proper rigid body column from the [Fx2] matrix given by 
equation ( 78 ) 


Ss/A 


1 for a symmetrical analysis. 

= -1 for an antisymmetrical analysis. 


Generalized Aerodynamic Force Coefficients 


The generalized response force coefficients are read from the Doublet Lattice program in 
the form of complex matrices calculated at each <0 value of interest. Relating these 
coefficients to the [M4J and [M5] matrices provides: 


[M4] = -2q [QR]j^^ 

[M5I - 


( 86 ) 


The generalized excitation forces are formed by again extracting the proper rigid body 
columns from the generalized response force coefficients. Since there is no gradual 
penetration, the matrix degenerates to a {C3} matrix. For <0 = 0, the generalized 
excitation is: 

{Qe} = {C3} ^ 


and for <0 > 0: 

{Qe} - {C3} = 28s/a ^ {Q'} 


( 88 ) 
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where: 


{Q'} = a complex column matrix of generalized aerodynamic force coefficients 

formed by extracting the proper rigid body column from [QR] matrix. 

4.9 COORDINATE TRANSFORMATION AND AERODYNAMIC AXIS SHIFTING 

All interpolation that is performed in L217 (EOM) using the SA arrays from L215 
(INTERP) is carried out in the local structural axis system of each aircraft component 
(sec. 4.7 of ref. 6). The aerodynamic geometry data is input in the reference axis system 
and, therefore, it must be transformed into the local structural axis for interpolation. 
This transformation is complicated by two factors. First, in most cases the origin of the 
local structural axis system and the origin of the local aerod 3 Tiamic axis system are not 
the same. Second, due to the requirements for good aerodynamic modeling, certain 
components may have to be placed outboard or inboard, above or below, or fore or aft, of 
their actual locations. Therefore, for proper interpolation, the reference axis system 
locations of points on such a surface must be corrected in order to reflect their true 
orientation. 

In order to account for the difference, the transformation equation for each component 
has coordinate shift values which can be used to define the correct position of any 
surface points. The transformation equation for any point on a surface is: 


X' 


r 

X' 


XoA 


■Xsh' 


y 

= [R] 


Y 

- ' 

Yoa 

+ ' 

Ysh 


.Z. 



Z. 


ZOA- 


■ZSH‘ 



^ s 

h 


RAS 


RAS 


RAS- 


where: 


y/s» ^/s 

XraS. YraS. Zras 

^OAraS’ ^OAras 

^SHrAS’’^SHras«ZsHras 

[R] 


= Local structural axis system location of any point 
on a surface 

= Reference axis system location of the point 

= Reference axis system location of the origin of the 
aerodynamic axes of the surface 

= Reference axis system shift values that will 
properly orient the aerodynamic axes to the 
structural axes 

= Rotation matrix that rotates the coordinates from 
the reference axis system into the local structural 
axis system 
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The rotation matrix is defined as: 



0 

0 

COS0 

sind 

-sin0 

cosd 


( 90 ) 


with 0 being the surface dihedral angle. 

The reference axis system coordinates of the point and the origin of the aerodynamic 
axis system, and the surface dihedral angle are obtained from the aerodynamic 
geometry data read by L217 (EOM). The shifting values are input by the user for each 
surface. 

To understand the calculation of the shifting values, consider the horizontal tail shown 
in figure 13. Here, a typical aft body modeling problem is illustrated. Both FLEXSTAB 
and Doublet Lattice require interference bodies to be of constant cross-section. Thus, in 
the aft portions of the fuselage where in reality the body tapers, the interference body 
must be kept constant. This could place the horizontal tail outboard, forward, aft, above 
or below its actual location. Just calculating the differences between the local 
aerodynamic and structural axis system origin would not reflect the true position of the 
aerodynamic points with respect to the interpolation axes (which are the local 
structural axes). Therefore, a set of correction values are needed. 

Taking into account these corrections, the equation for the shift values is: 


'Xsh' 


'Xoa' 

YSH 

= ^ 

Yoa 

.ZsH. 


.ZoA. 


RAS 

where: 



'Xos 


^Xc 

- 

YOS 

+ ■ 

Yc 

RAS 

.Zos. 

RAS 

zc. 


(91) 


XoStjac’^OS ac»^OSt>a = Reference axis system location of the structural 

rt-AiO x\.A.o fCzVo j 1 n ji Cl A 

axis system used for the SA array. 

^^RAS’ ^^RAS’ ^^RAS -Desired correction values to apply to the 

aerodynamic points. These values are defined in 
the reference axis system and are positive in the 
positive reference axis system directions measuring 
from the aerodynamic axis system to the structural 
axis system. 

The proper calculation of the shifting values requires the user to know the location of 
the local structural axis system origin used in L215 (INTERP), and the local 
aerodynamic axis system origin established by FLEXSTAB or Doublet Lattice. This last 
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Figure 13. — Aerodynamic Axis Shifting 


value for FLEXSTAB aerodynamics is defined by the analyst when setting up the 
FLEXSTAB idealization. Doublet Lattice, however, defines the aerodynamic origin 
internally. In Doublet Lattice, the aerodynamic origins are defined in the following 
way: 

• For slender bodies, the origin is the leading edge location of the first segment. 

• For thin bodies (or panels), the origin is the XI, Yl, Z1 coordinates of the panel. 
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5.0 PROGRAM STRUCTURE AND DESCRIPTION 


Program L217 (EOM) is structured as an overlay system consisting of a main overlay 
and four primary overlays listed below and displayed in figure 14. 

(L217,0,0) L217VC 

(L217,2,0) STRUCT 

(L217,3,0) FLXAIC 

(L217.4.0) DLAIC 

(L217,5,0) DLPRES 


5.1 PROGRAM L217VC 

The main overlay ,(L217, 0,0), is named L217vc (the "v” and "c” are version and 
correction identifiers). It reads data cards that direct the execution of the primary 
overlays and also aids communication between the primary overlays via labelled 
common blocks. 


5.2 PROGRAM STRUCT 


The 2,0 primary overlay, STRUCT (L217,2,0), generates the equations of motion 
matrices that are dependent upon the structural data only. 

[Mx] = Generalized stiffness 

[M 2 ] = Generalized damping (optional) 

[M 3 ] = Generalized inertia 

Depending on the options chosen, STRUCT may read from several magnetic files: 
GSTIFF A generalized stiffness matrix [K]. ([K] = [Ml],. ) 

■*U 3. 

GMASS A generalized mass matrix [M]. ([M] = [M 3 ] ) 

Ua 

SATAP Arrays of modal data to be used in the generation of control surface freedoms to 

be added to [K] and [M]. 

Because certain constants (q and Vt) and defined in STRUCT, and because the 
matrices generated by STRUCT are required in each of the aerodynamic overlays, 
STRUCT must be executed prior to any of the aerodynamic overlays. 
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Generates structural Uses F LEXSTAB 

equations of motion AlC matrices 

matrices [M J , [Mo] , to: 

and [M^l 


Uses Doublet 
Lattice AlC 
matrices to: 


Uses Doublet Lattice 
generalized forces 
and unsteady pressures 
to: 


Generate aerodynamic equations of motion matrices 
([M^l , [Mg] , and [^] ), plus the matrices of 
aerodynamic forces per surface ([F-j] , IF2], and [FG] ). 


! 

i 


I Figure 14, - L217 (EOM) Overlay Structure 
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5.3 AERODYNAMIC DATA PROCESSING 


The remaining primary overlays generate the equations of motion matrices that are 
dependent upon the aerodynamic data. The matrices are: 

[FREQM] = Array of frequencies at which the aerodynamics are defined 


[M4] 

[M5] 


= Generalized aerodynamic stiffness matrix 
= Generalized aerodynamic damping matrix 

= Array distance(s) from the gust reference point to the gust control 
point(s) of the gust zones 


= Gust forcing function matrix 

The matrices are written onto the file *'EOMTAP” following the structural matrices 
generated by STRUCT. 


In addition, the aerodynamic overlays produce the aerodynamic force matrices per 
surface and write them onto the file "EOMLOD” for subsequent use in L218 (LOADS) 
(ref. 9). The file includes the following matrices for each surface: 

= Local axis coordinates of surface node points 

= Area of surface elements 

= Real part of aerodynamic response forces 

= Imaginary part of aerodynamic response forces 

= Gust forces 

The last three matrices are repeated for each frequency at which they are defined. 


[x.y] 

{A} 

[Fi] 

[F2] 

[FG] 


5.3.1 PROGRAM FLXAIC 

The 3,0 primary overlay, FLXAIC (L217,3,0), generates the aerodynamic equations of 
motion and force matrices using the FLEXSTAB (ref. 7) aerodynamic influence 
coefficient (AIC) matrix. FLXAIC reads from the files: 

"GDTAPE” FLEXSTAB geometry data. 

"SDINDX” Table describing the file of aerodynamic data, SDDATA. 

"SDDATA” AIC matrix [Aj^ for a symmetric or antisymmetric case. 

"SATAP” Arrays of modal data. 


SO 



5.3.2 PROGRAM DLAIC 


The 4,0 primary overlay, DLAIC (L217,4,0), generates the aerodynamic equations of 
motion and force matrices using the Doublet Lattice (ref. 8) AIC matrices. DLAIC reads 
from the files: 

”NGETP” Doublet Lattice geometry data. 

"*NAETP” Doublet Lattice aerodynamic data ([F], [D] etc.). 

"SATAP” Arrays of modal data. 

5.3.3 PROGRAM DLPRES 

The 5,0 primary overlay, DLPRES (L217,5,0), generates the aerodynamic equations of 
motion and force matrices using the Doublet Lattice (ref. 8) generalized forces and 
unsteady aerodynamics. DLPRES reads from the files: 

”NGETP” Doublet Lattice geometry data. 

"NAETP” Doublet Lattice aerodynamic data ([Q],[ACP], etc.). 

Although L217 (EOM) serves as a module of the DYLOFLEX system, it can be operated 
as a standalone program. L217 (EOM) is heavily dependent upon data generated prior 
to its execution. Figure 15 displays the different files that the program may require, 
according to user specified options. When the program is run by itself, it becomes the 
user’s responsibility to generate the files of input data in the format required by the 
program. 
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FLEXSTAB 



GD1APE NGETP GSTIFF GMASS SATAP 



To L221 (TEV156) To L218 (LOADS) 

or L219 (EQMOD) 


Figure 15. — L2 17 (EOM) Communication Via External Files 







6.0 COMPUTER PROGRAM USAGE 


The program was designed for use on the CDC 6600. The machine requirements to execute 
L217 (EOM) are: 


Card reader Read control cards and card input data. 

Printer Print standard output information, optional intermediate calculations, 

and diagnostic messages. 

Disk storage All magnetic files not specifically defined as magnetic tapes are 

assumed to be disk files. 

Tape drive For permanent storage of data. Magnetic files are copied to and from 

magnetic tapes with control cards before and after program execution. 


The program L217 (EOM) is written in FORTRAN and may be compiled with either the RUN 
or FTN compiler. L217 (EOM) may be executed on either the KRONOS 2.1 or NOS operating 
system. 

6.1 CONTROL CARDS 


The following list is a typical set of control cards used to execute L217 (EOM) using the 
absolute binaries from the program’s master tape. 

Job Card 
Account Card 




• 

REQUEST(MASTER,F=I,LB=KL,VSN=66XXXX) 

REWIND(MASTER) 

SKIPF(MASTER) 

COPYBF(MASTER,L217) 

RETURN(MASTER) 


L217. 

• 

EXIT. 

DMP(0,field length) 
— End-of-record 


"Retrieve the 
program from 
.its master tape 


{ Prepare input 
data files 

{Execute L217 (EOM) 

f Save output 
data files 


Card Input 
Data 


— End-of-file 
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The following list is a typical set of control cards used to execute L217 (EOM) using the 
relocatable binaries from the program’s master tape. 


Job Card 
Account Card 

• 

e 

'Retrieve the 
program from 
its master 
tape 

{ Prepare input 
data files. 

( Retrieve DYLIB, 
the DYLOFLEX 
alternate library 

LDSET(LIB=DYLIB,PRESET=INDEF) {Load L217 (EOM) 

LOAD(REL217) 

NOGO. 

RETURN(REL2 1 7 ,D YLIB) 

L217. {Execute L217 

• r Save output 

• [data files 

• 

EXIT. 

DMP(0, field length) 

— End-of-record 


# 

REQUEST(MASTER,F= I,LB= KL,VSN= 66XXXX) 
REWIND(MASTER) 

SKIPF(MASTER,2) 

COPYBF(MASTER,REL217) 

RETURN(MASTER) 

• 


Card Input 
Data 


— End-of-file 
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6.2 RESOURCE ESTIMATES 


The computer resources used (core requirements, tapes, printed output, central 
processor seconds, etc.) are a function of the problem size and the program options used. 

Field Length 

The field length (core) required by L217 (EOM) is dependent upon the problem size and 
the program module(s) used. Core must be requested for the largest module run, either 
the structural data overlay or one of the aerodynamic overlays. The structural and 
aerodynamic overlays (primaries) are broken into even smaller chunks called secondary 
overlays. Therefore, core must be requested for the longest secondary overlay to be 
executed. The equation to calculate the field length of an overlay is 

Field Length = FLp + FL^ 

where FLp is the program length, and FL^ is the array storage required. 

The following pages contain FLp and a simplified formula for FL^ for all L217 (EOM) 
secondary overlays. The user must determine which ones will be executed, which is the 
longest for the problem to be run, and the required field length. 

It should be noted that FLp is given as an octal base number. Normally, FLa will be 
calculated as a decimal number and converted to octal before being added to FLp. 

The structural overlay is relatively short. One of the aerodynamic overlays will usually 
be the largest. 


The following abbreviations are used in the field length formulas: 
NDOF( 2 ;s Number of control surface freedoms being added (card 2.1). 
Number of original structural freedoms (card 2.1). 


NDOFgT 

NMODES 

Mmax 

NELE 

NGCP 

NIBFP 

NSBFP 


Number of modes defined (number of generalized coordinates total) = 

NDOFgx + NDOF^g. 

Minimum of 20 or NMODES. 

Number of control surface masses and inertias added to the structural 
matrices (card 2.9.1). 

Number of gust control points (card 3.9.1). 

Number of interference body force points. 

Number of slender body force points. 



NTBFP 


Number of thin body force points. 


NAFP Number of output aerodynamic force points = NTBFP + NSBFP. 

NBOX Number of aerodynamic boxes = NTBFP + NIBFP. 

NTOTF Number of aerodynamic elements total = NTBFP -h NIBFP -h NSBFP. 

Nmax Maximum number of nodes on a surface. 

Primary Overlay L217,2,0 
Program Name STRUCT 


If control surface freedoms are added, FORMAC will be longest. 
If not, FORM13 will be the longest. 


Secondary 

Overlay 

Program 

Name 

FLp 

(octal) 

FLa 

(decimal) 

L217,2,l 

RSTRCT 

35000 

6000 

L217,2,2 

FORMAC 

74000 

NELE* (NMODES + 1) + 600 
+(NDOFcs*NDOFst) 

L217,2,3 

FORM13 

47000 

NMODES^ + NMODES 


Primary Overlay L217,3,0 
Program Name FLXAIC 


For large problems L217,3,4-FLXFRC will be the longest. 
For small problems L217,3,3-FLXMOD will be the longest. 


Secondary 

Overlay 

Program 

Name 

FLp 

(octal) 

FLa 

(decimal) 

L217,3,l 

FLXCDS 

36000 

0 

L217,3,2 

FLXGEO 

66000 

0 

L217,3,3 

FLXMOD 

73000 

3*NAFP*(Mmax + 2) 

L217,3,4 

FLXFRC 

42000 

4*NAFP*(Mmax + D + 5* NTOTF 

L217,3,5 

FLXODR 

42000 

NAFP*(Mmax + 5) + NMODES*Nmax 

L217,3,6 

FLXEOM 

42000 

2* NMODES* (NMODES +NGCP) 
+ NAFP*(Mmax + D 
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Primary Overlay L217,4,0 


Program Name DLAIC 

For large problems, L217,4,6 -SOLVDW is usually the longest. 
For small problems, L217,4,3 -DLAMOD is usually the longest. 


Secondary 

Overlay 

Program 

Name 

FLp 

(octal) 

FLa 

(decimal) 

L217,4,l 

DLACDS 

37000 

0 

L217,4,2 

DLAGEO 

65000 

0 

L217,4,3 

DLAMOD 

75000 

Mmax*(3*NAFP+NTBFP)+5*NAFP 

L217,4,4 

CPANDW 

40000 

NTBFP»(Mmax+2*NMODES) 
+ 2*(NAFP+NMODES) 

L217,4,5 

FORMWS 

41000 

2*NBOX*NMODES 

L217,4,6 

SOLVDW 

52000 

NBOX*(NBOX+4*(NMODES+NGCP)+3)^ 

L217,4,7 

DLAFRC 

41000 

2*((NSBFP*NMODES)+NBOX 
+ (NAFP*Mmax» +5*NT0TF 

L217,4,8 

DLAODR 

43000 

2*(NMODES*Nmax+NAFP*Mmax) 

+6*NAFP 

L217,4,9 

DLAEOM 

43000 

2*NMODES*(NMODES+NGCP) 

+NAFP*(3’*‘Mmax+l) 


^The formula given above for SOLVDW is to perform a solution in core. SOLVDW will 
partition the problem and run with as little as FLa = NBOX^CNMODES+NGCP) 
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Primary Overlay L217,5,0 
Program Name DLPRES 


DLPEOM will normally be the longest. 


Secondary 

Overlay 

Program 

Name 

FLp 

(octal) 

FLa 

(decimal) 

L217,5,l 

DLPCDS 

35000 

0 

L217,5,2 

DLPGEO 

63000 

0 

L217,5,4 

DLPFRC 

40000 

7*NTOTF + 3*NAFP 

L217,5,5 

DLPODR 

42000 

Nmax* (NMODES + 4) + 8* N AFP 

L217,5,6 

DLPEOM 

41000 

4*NMODES^ + 2*NMODES 


Time Estimate 

The central processing time required (CP seconds) is dependent upon the options chosen 
and the problem size. The total time will be the sum of the time required by each 
overlay executed: 

+ CPpLXAIC + CPdlAIC + CPdlPRES 
NFREQM 

CPdLPRES ^ 5 * NFREQM 
where: 

NFREQM = Number of frequencies at which the aerodynamics are defined 
NBOX = Number of boxes on the thin and interference bodies 

Printed Output 

The amount of printed output is dependent upon card options (cards 2.4, 2.6, 2.7 and 
3.11), the size of the arrays printed, and the number of frequencies at which 
aerodynamic data is processed. At compilation time, the print line limit is set to 40 000 
which should be sufficient for any single L217 (EOM) data case. 


CPtoT 

CPsTRUCT 

cpflxaic 

cPdlaic 


= CPstRUCT 
« 10 
10 


(l 


. , 520X ‘1. 


\ 40 
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6.3 CARD INPUT DATA 

The task(s) performed by L217 (EOM) are broken into several subtasks, each with its 
own section of code known as a primary overlay. The entire set of primary overlays is 
driven by a small program (main overlay) named L217vc. 

L217vc reads program directive cards to: 

• Assure that the data being read is intended for EOM 

• Determine which section of code (primary overlay) of EOM is to be executed next. 

Each primary overlay reads one set of card input data. A listing of the primary overlays 
and the card set numbers follows. The same card set (3) is read by all overlays 
processing aerodynamic data. 


Overlay 

Name 

Purpose of Primary Overlay 

Card Set 
Read 

STRUCT 

Generate [Mj], [M2] equations of motion 

matrices based upon structural data. 

2 

FLXAIC 

Generates the aerodynamic coefficient matrices for the 
equations of motion based on the quasi-steady 
aerodynamic assumption and using the aerodynamic 
data from FLEXSTAB (ref. 7) 

3 

DLAIC 

Generates the aerodynamic coefficient matrices for the 
equations of motion based using Doublet Lattice AIC^s 
(ref. 8). Either quasi-steady or full unsteady 
aerodynamics can be used. 

3 

DLPRES 

Generates the aerodynamic coefficient matrices for the 
equations of motion using the Doublet Lattice 
generalized forces (ref. 8) for full unsteady 
aerodynamics. 

3 


The order in which these card sets are input (and the overlays are executed) is 
displayed in figure 16. The overlay STRUCT must be executed before any of the 
aerod 3 mamic overlays. 
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Format of Card Input Data 

All card data is read in fixed fields (specific columns of the cards). On the pages that 
follow, the required card columns are defined next to each keyword or variable. The 
following conventions are used throughout the program: 

• All floating point variables are read with format ElO.O. 

• All integer variables are read with the format 15. 

• All hollerith variables (keywords, etc), are read with the format AlO. 

Therefore, all data fields end on a card column that is a multiple of five. 

When the program is trying to recognize keywords, it checks only the first four 
characters. 

6.3.1 PROGRAM L217 

The driving module of L217 (EOM) (the main overlay itself is also known as L217vc) 
will read and check the first four characters of the input card to determine the type of 
execution that is to take place. 

Note: All underlined capital characters contained in the KEY WORD/VARI ABLE field 
of the input card sets must be left-justified and punched in the card columns 
specified in the COLS, field of the input card sets. 


L217 (EOM) Program Directive Cards 

All program directive cards have a in column one, the first character of the four 
character keywords. 

Card Set. 1.0-Introduce L217 (EOM) Data 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

$EOM 

AlO 

Keyword introducing EOM data. This must be the first 
card read by EOM. 

11-70 



Available for comments 
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Card 1.1-Title Card Labelling Printout 

Optional card to be omitted or repeated as often as desired before any program 
directive card except $EOM. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

$TIT Te 

AlO 

Keyword indicating a title card to be read, printed, 
and ignored. 

11-70 


1 

1 

Available for cominents 


Card 1,2-Indicate a Program Checkout Run 

Optional card used only for program checkout purposes. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

$NON-stop 

AlO 

Keyword telling the program to try continuing execution 




even if fatal errors are diagnosed. 

11-70 



Available for comments 


6.3.2 STRUCTURAL CARD INPUT DATA 

After reading card 1.0, (cards 1.1 and 1,2 are optional), the next card that must be read 
is card 2.0, which introduces the structural input data. Cards 2.1 through 2.3 define the 
problem size, the type of analysis, and the type of generalized mass and stiffness data 
input. Card 2.4 is optional and is used to print out intermediate data for problem 
checkout purposes. Card 2.5, defining the modal input, is only required if control surface 
inertial coupling is to be calculated. Cards 2.6 through 2.8 define the generalized 
inertial, stiffness, and damping (if calculated in L2 1 7 (EOM)). Finally, card 2.9 defines 
the control surface inertial and nodal data needed for the inertial coupling calculations. 

Once the processing of the structural data has been completed, the program reads 
another card to determine what is to be done next. The user may specify an 
aerodynamic data set, or new structural data ($STR - card set 2), or may terminate the 
program’s execution with card set 4 ($QUIT). 

All the cards in the set are optional except for cards 2.0 ($STR), 2.1 (SIZE), and 
2.10 (END). Card 2.1 must follow card 2.0 and card 2.10 must terminate the data set 
and precede the next program directive card. The flow of cards is displayed in figure 17. 

For any optional card (2.2 through 2.9) that is omitted, the keywords or variables 
appearing on the card will assume their default values. 
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Card Set 2.0-Introduce Structural Card Input Data 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

$STRuctural 

AlO 

Keyword introducing the structural card input data 
and causing the STRUCT overlay to be called into 
execution* 

The structural card input data (Cards 2.1-2.10) must 
appear between $STR and the next program directive 
card (either $TIT, $CHE, or $AER) . 


11-70 



Available for comments 


Default values have been established for most of the items input in the structural data 
set. All cards in the set are optional, with the exception of 2.0 ($STR), 2.1 (SIZE), and 
2.10 (END). Card 2.1 must follow 2.0, and 2.10 must terminate the structural data and 
precede the next program directive card. Cards 2.2 through 2.9 may be in any order. 

Card 2.1-Size of Equations of Motion Matrices 

Card 2.1 must follow card 2.0. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SIZE 

AlO 

Keyword indicating a "SIZE" card. 

11-15 

NDOFST 

15 

The number of structural degrees of freedom. 

The size of the input generalized mass and stiffness 
matrices. 

16-20 

NDOFCS 

15 

The number of control surface degrees of freedom to be 
added to the structural freedoms. 


NOTE: The equations of motion matrices generated by L217 (EOM) will contain NDOF 
degrees of freedom where NDOF = NDOFST + NDOFCS. If the control surface 
freemdoms were accounted for in the generation of the structural data, then 
NDOFCS = 0 and NDOFST is set equal to the total number of freedoms coming 
from the vibration program. 
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Card 2.2-Symmetric/ Antisymmetric Analysis 


Card 2.2 is optional. The default analysis type is symmetric. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

ISYM 

AlO 

One of two keywords may be used to choose the type of 
analysis to be performed by L217 (EOM) . 

SYMMetric 
ANT I - symme tr ic 
Default: SYMM 

11-70 



Available for comments 


Card 2.3-Full/Half Airplane 

Card 2.3 is optional. By default, the input generalized mass and stiffness matrices 
are assumed to be for a full airplane. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

o 

1— i 
1— 1 

IHALF 

AlO 

One of two keywords may be used to specify that the input 
generalized mass and stiffness are for full or half airplane. 
FULL 
HALF 

Default: FULL 

11-70 



Available for comments 



Card 2.4-Checkout Printout 


Card 2.4 is optional. 



DESCRIPTION 


Keyword requesting the printing of intermediate 
data calculated during the structural data processing. 


Option indicating the amount of intermediate data to 
be printed. 

=1 Print message when each overlay is entered plus 
[ 4 )] , interpolated mode shapes at control surface 

nodes, [MCS] , and the final [Mi] ,[M 2 ] , and [M 3 ] . 

=3 Print all items listed under 1 plus control surface 
data summary, [A] , and [C] 

=5 Print all items listed under 1 and 3 plus [(j)^,^] , 

OX 

[(|) ] , and a message each time a VARDIM routine 

v^o 

is called. 


Card 2.5-File Containing the Input Mode Shapes 

Card 2.5 is optional. The default file name is ''SATAP. 


Keyword indicating a card defining the modal input data 
file name. 


Name of the file (tape or disk) containing the modal 
input data which will be used to calculate the inertia 
coupling effects for the added control surface freedoms. 

The name must start in column 11, begin with a letter, and 
have less than seven characters. 


Default: SATAP= " SATAP " 
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Card 2.6-Generalized Stiffness Matrix Source 


Card 2.6 is optional. Default values are shown below. 



Keyword indicating stiffness matrix data 


Name of the file containing the generalized 
stiffness matrix [K ] , Either "CARDS" or a 
file name starting in column 11, beginning with 
a letter, and having less than seven characters. 
If KTAPE="CARDS" the matrix will be read from 
cards (see Card 2.6.1) and the next two variables 
are meaningless. 

Default: KTAPE="GSTIFF" 


Number of logical files to skip on KTAPE before 
reading [K] - 
Default: KFILE=0 


Number of matrices to skip on KTAPE before 
reading [K] . 

Default: KMAT=0 


Option requesting the printing of [K] as read. 
Default: Do not print except when KTAPE=CARDS 


Card 2.6.1-Rows of Generalized Stiffness Matrix on Cards 
Card 2.6.1 is input only if KTAPE = CARDS on card 2.6 
Repeat card 2.6.1 for each row of [k1; i = 1, NDOFST. 


KEYWORD/ 

COLS . VARIABLE FORMAT 


1-10 I K(i,l) 



11-20 K(i,2) 


DESCRIPTION 


7E10.0 The (i,l) element of [K] the generalized stiffness 


The (i,2) element of [K] 



K(i, j) 


K(i, NDOFST) 



The (i,j) element of [K] 


Repeat the card with 7 numbers per card until all 
elements of row i are defined. 


































Card 2.7-Generalized Mass Matrix Source 


Card 2.7 is optional. Default values are shown below. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10. 

MASS 

AlO 

Keyword indicating mass matrix data. 

11-20 

MTAPE 

AlO 

Name of the file containing the generalized 
mass matrix fM] . Either "CARDS” or a file 
name - starting in column 11, beginning with 
a letter, and having less than seven characters. 
If MTAPE="CARDS" the matrix will be read from 
cards (see Card 2.7.1) and the next two variables 
are meaningless. 

Default: MTAPE=”GMASS" 

21-25 

MFILE 

15 

Number of logical files to skip on MTAPE before 
reading [M] . 

Default: MFILE=0 

26-30 

MMAT 

15 

Number of matrices to skip on MTAPE before reading 
[M] . 

Default: MMAT=0 

31-40 

PRINt 

AlO 

Option requesting the printing of [MJ as read. 
Default: Do not print except when MTAPE=CARDS 

41-70 



Available for comments 
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Card 2.7.1-Rows of Generalized Mass Matrix of Cards 


Card 2.7.1 is input only if MTAPE = CARDS in card 2.7. 
Repeat card 2.7.1 for each row of [M], i = 1, NDOFST. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 


7E10.0 

The (i,l) element of [M] , the generalized mass matrix. 

11-20 

M(i,2) 


The (i,2) element of [M] 

• 



The (if j) element of [M] 


M(i, NDOFST) 


Repeat the card with 7 numbers per card until all 




elements of row i are defined. 


Card 2.8-Introduce Structural Damping 

Card 2.8 is optional. The default is no viscous structural damping. 


COLS. 

REWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

DAMPing 

AlO 

Keyword introducing the structural damping factors 
following in Card 2.8.1 




Available for comments 


Card 2,8,1-Structural Damping Factors 

Card 2.8.1 must follow card 2.8. It is repeated until all NDOFST terms are defined. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

gSD (1) 

7E10.0 

Element 1 of {gSD} 

11-20 

gSD (2) 


Element 2 of {gSD} 


gSD (i) 


Element i of { gSD} 


gSD (NDOFST) 


Repeat the card with 7 numbers per card until 




all elements are defined. 
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Card 2.9-Introduce the Addition of Control Surface Freedoms 


Omit cards 2.9 and 2.9.1 if NDOFCS = 0 on card 2.1, 

Repeat cards 2.9 and 2.9.1 for each control surface affected by the additional 
freedoms. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SQRFace 

AlO 

Keyword introducing a new surface. 

11-15 

IS 


Modal Interpolation Program surface number 
(Reference 5) specifying where mode shapes for 
this surface are to be found. 

16-20 

NNODEI 

15 

Number of nodes (mass points) on this control 
surface. 

21-30 

Trans 

AlO 

A "T” in column 21 requests a translational 
freedom for each mass point. 

31-40 

Rotl 

AlO 

An "R“ in column 31 requests a streamwise rotational 
freedom, e , for each mass point, (local aerodynamic axes) 

41-50 

Rot2 

AlO 

An "R" in column 41 requests a rotational freedom, 
, which is perpendicular to the freestream, for 
each mass point, (local aerodynamic axes) 


Local 

AlO 

An "L" in column 51 indicates that the coordinates on 
Card 2.9.1 are in the local structural axis system. 
Default: Reference Axis System 


Card 2.9.1-Mass Data for Surface IS 
Card 2.9.1 must follow card 2.9. 

It is repeated NNODEI times, once for each node on the control surface. 


Hi 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

Xi 


The coordinates of node i on surface IS. The 
coordinates are in the reference axis system 
unless specified as Local on Card 2.9. 

11-20 

Yi 

QBQDH 

21-30 

Zi 


31-40 

MASSi 


The mass at node i. 
Whole airplane. 

41-50 

INERTYi 

I^QQQII 

The mass moments of inertia node i about the local aerodynamic 
y and x axes. Whole airplane. 

Note: INERTYi is required only if R0T1=R on Card 2.9. 

INERTXi is required only if ROT2=R on Card 2.9. 

51-60 

INERTXi 


61-70 



Available for comments 
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Card 2.10-End Structural Card Input Data 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

END 

AlO 

Keyword indicating the end of structural card 
input data. 

11-70 



Available for comments 


6.3.3 AERODYNAMIC CARD INPUT DATA 

After card set 2 (2.0 through 2.10) and the processing of structural data, the next card 
read by L217 (EOM) must introduce aerodynamic card data. Card 3.0 introduces the 
aerod 3 mamic card data and indicates the source of the AIC’s or generalized forces, 
FLEXSTAB or Doublet Lattice. Some of the cards within card set 3.0 (3.13 through 
3.16) vary according to the source of the aerodynamic data used. 

Once the processing of the aerodynamic data has been completed, the program reads 
another card to determine what is to be done next. The user may specify another 
aerodynamic data set, or new structural data ($STR-card set 2), or may terminate the 
program’s execution with card set 4.0 ($QUIT). 

Default values have been established for many of the items defined in the aerodynamic 
card set. All cards in the set are optional with the exception of3.0($AER), 
3.1 (CONSTANTS), and 3.17 (END). Card 3.1 must follow 3.0 and card 3.17 must 
terminate the data set and precede the next program directive card. Optional cards 3.2 
through 3.8 need only be present if the analyst wishes to use file names and/or spacing 
parameters that are different from the default values appearing in the data card 
write-up. The flow of cards is displayed in figure 18. 




Required cards Optional cards 



Figure 18. — L217 (EOM) Aerodynamic Card input Data, Card Set 3.0 






KEYWORD/ 

VARIABLE 


FORMAT 


DESCRIPTION 


COLS. 


1-10 

$AER 

— 

AlO 

Keyword introducing the set of cards describing 
the aerodynamic data processing. Check the next 
two words for aerodynamic type. 

11-20 

ISOURC 

AlO 

Keyword indicating the source of the aero data. 
There are two possibilities: 

FLEXstab - FLEXSTAB 
DOUBlet - Doublet Lattice 

21-30 

ITYPE 

AlO 

Keyword indicating the type of aero data to be used. 
There are two possibilities: 

AIC - AIC matrices 

PRESsure - Generalized forces and pressures 
Note: If the source is FLEXSTAB the type must be 

AIC . The PRES option is valid only with Doublet Lattice. 

Lattice. 

31-70 


J 

Available for comments 




















Card 3.1-Aerodynaniic Constants 


Card 3.1 must follow card 3,0. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

CONStants 

AlO 

Keyword indicating a card of aerodynamic constants. 

11-20 

RHO 

ElO.O 

p - density of air at altitude. 

21-30 

VEL 

ElO.O 

= velocity (true airspeed) [units of length per second] 
Note: q, the dynamic pressure, is calculated as 

q = *SPV‘ 

31-35 

NTDEF 

15 

Indicates the number of the mode corresponding to 
the translation (vertical or sideO freedom. That 
mode will be used to find the gust forces when 
using Doublet Lattice pressures at a non-zero 
reduced frequency. 

36-40 

NPDEF 

15 

Indicates the number of the mode corresponding to the 
pitch or yaw freedom. That mode will be used to find 
the gust forces when using Doublet Lattice pressures at 
a reduced frequency of zero. 


NTDEF and NPDEF are not used, except with Doublet Lattice generalized forces 
and pressures. 

Card 3.2-File Containing the Input Mode Shapes 

Card 3.2 is optional. The default file name is ''SATAP.” 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SATAp 

AlO 

Keyword introducing replacement name for the modal 
interpolation data file. 

11-20 

ISATAP 

AlO 

Name of the file (tape or di4k) containing the modal 
input data. 

The name must start in column eleven, begin with a 
letter, and have less than seven characters. 

Default: ISATAP = "SATAP" 

21-70 



Available for comments 
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Cards 3.3 and 3.4-FLEXSTAB Input Data Files 


Omit cards 3.3 and 3.4 if Doublet Lattice aerodynamics are being used. The cards 
are optional even if FLEXSTAB aerodynamics are being used. Default file names 
are shown below. 


Card 3.3-File Containing FLEXSTAB Geometry Data 

Card 3.3 is optional. The default file name and case number are "GDTAPE” and 1, 
respectively. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

GDTApe 

AlO 

Keyword introducing replacement name for the 
FLEXSTAB geometry data file. 

11-20 

GDTAPE 

AlO 

Name of the file (tape or disk) containing the 
FLEXSTAB geometry data. 

The name must start in column eleven, begin with 
a letter, and have less than seven characters. 

21-25 

I GCAS E 

15 

FLEXSTAB geometry data case number from which L217 
(EOM) must extract geometry information. 

Default: IGCASE= 1 

26-70 



Available for comments 
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Card 3.4-Files Containing FLEXSTAB Aerodynamic Data 


Card 3.4 is optional. The default file names are ”SDINDX” and "SDDATA.” 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SDSStp 

AlO 

Keyword introducing replacement name for the FLEXSTAB 
aerodynamic data file. 

11-20 

SDINDX 

AlO 

Name of the file (tape or disk) containing the index 
table to the FLEXSTAB aerodynamic data. 

The name must start in column eleven, begin with a 
letter, and have less than seven characters. 

Default: SDINDX = "SDINDX" 

21-30 

SDDATA 

AlO 

Name of the file (tape or disk) containing the FLEXSTAB 
aerodynamic data. Note: SDINDX SDDATA. 

The name must begin in column 21, begin with a letter, 
and have less than seven characters. 

Default: SDDATA = "SDDATA” 

31-35 

ISCASE 

15 

FLEXSTAB aerodynamic data case number from which L217 (EOM) 
must extract [A-. ] 

Default: ISCASE = 1 

36-70 



Available for comments 
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Cards 3.5 and 3.6-Files Containing Doublet Lattice Data 

Omit cards 3.5 and 3.6 if FLEXSTAB aerodynamics are being used. 

The cards are optional even if Doublet Lattice aerodynamic data is being used. 
Default file names are shown below. 


Card 3.5-File Containing Doublet Lattice Geometry Data 

Card 3.5 is optional. The default file name is ”NGETP” and the default 
IGCASE =1. 


COLS. 

KEYWORD/ 

VARIABLE 

FORIvIAT 

DESCRIPTION 

1-10 

NGETp 

AlO 

Keyword introducing a replacement name for the 
input Doublet Lattice geometry file. 

11-20 

NGETP 

AlO 

Name of the file (tape or disk) which contains 
the Doublet Lattice geometry data. 

21-25 

IGCASE 

15 

Number of the data case (logical file) on "NGETP" 
to be used. 

26-70 



Available for comments 


Card 3.6-File Containing Doublet Lattice Aerodynamic Data 

Card 3.6 is optional. The default file name is ”NAETP” and the default case 
number is 1. 


_COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

NAETp 

AlO 

Keyword introducing a replacement name for the 
input Doublet Lattice Aerodynamic file. 

11-20 

NAETP 

AlO 

Name of the file (tape or disk) which contains the 
Doublet Lattice aerodynamic data. 

21-25 

lACASE 

15 

Number of the data case (logical file) on "NAETP" 
to be used. 

26-70 



Available for comments 
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Card 3.7-File to Contain the Equations of Motion Coefficients 


Card 3.7 is optional. The default name is "EOMTAP.” 


COLS. 

KETiTWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

EOMTap 

AlO 

Keyword introducing replacement name for the output 
equations of motion file. 

11-20 

EOMTAP 

AlO 

Name of the file (tape of disk) to contain the 
equations of motion matrices. 

The name must begin in column eleven, begin with a 
letter, and have less than seven characters. 

21-70 



Available for comments 


Card 3.8-File to Contain Aerodynamic Force Matrices 

Card 3.8 is optional. The default file name is "EOMLOD.’' 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

EOMLod 

AlO 

Keyword introducing replacement name for the output 
file to contain the force matrices required by L218, 
the Load Equations program. 

11-20 

EOMLOD 

AlO 

Name of the file (tape or disk) to contain the force 
matrices required by L218, the Load Equations program. 
The name must begin in column eleven, begin with a 
letter, and have less than seven characters. 

21-70 



Available for comments 


Card 3. 9- Doublet Lattice Quasi-Steady Option 

Card 3,9 is optional. The quasi-steady solution may be requested only with Doublet 
Lattice AIC aerodynamic data for a reduced frequency (k-value) of zero. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

QUASI 

AlO 

Keyword requesting a Quasi-Steady solution for the 
Doublet Lattice AIC option. 

11-70 



Available for comments 
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The quasi-steady option limits processing to the first reduced frequency solution 
from Doublet Lattice files. Therefore, the first reduced frequency must be zero. 


Card 3.10-FLEXSTAB Thin Body Control Point Option 

Card 3.10 is optional. It may be used only with FLEXSTAB aerodynamics. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 


1-10 

CENTroids 

AlO 

Keyword requesting L217 (EOM) to use the centroids 
of thin bodies' panels for both the force points and 
the control points. 

Default: Centroids are used for force points and control 

points will remain as defined in FLEXSTAB. 

11-70 



Available for comments 



Cards 3.11 and 3. 11.1 -Gradual Penetration 

Cards 3.11 and 3.11.1 are optional. They are used only when gradual penetration is 
desired. Gradual penetration can not be requested when using Doublet Lattice 
generalized forces and pressures. 

Card 3.11-Introduce Gradual Penetration Control Points 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

GRAD. pen. 

AlO 

Keyword introducing gradual penetration card input 
data. 

11-20 

XG 

EIO.O 

The gust reference point in the reference axis system. 

21-25 

NGCP 

15 

The number of gust control points. 
(liNGCP<35) 

26-70 



Available for comments 
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Card 3.11.1-Gradual Penetration Control Points 

Card 3.11.1 is input only after card 3.11. Card 3.11.1 must be repeated for each 
gradual penetration control point; i = 1,NGCP (from card 3.11). 


— 

COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

XGCP^ 

ElO.O 

The reference axis x-coordinate of the gust control 
point for zone i 

XGCP^_^^ > XGCP^ 

11-20 

XGLE . 
1 

ElO.O 

The reference axis x-coordinate of the leading edge 
of zone i. 

XGCP^ > XGLE^ , XGLE^_^j^ > XGLE^ 

21-70 



Available for comments 


The program will group all aerodynamic elements into the zones generated by the 
XGLE] of card 3.11.1. The gust will reach each panel of zone i at the time it 
reaches the gust control point XGCPi- XGLEj must be forward of the first 
element’s control point. 

Card 3. 12*- Request Printed Output 

Card 3.12 is optional. If omitted the following items are printed by default: Card 
input data and interpretation of card input data. Summary of geometry data read 
from GDTAPE or NGETP, 

Card 3.12 must be input once for each optional item to be printed. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

PRINt 

AlO 

Keyword introducing optional item to be printed. 

11-20 

ITEM 

AlO 

Name of the item to be printed chosen from the 
following list: 

MOD^S: Modal displacements at aerodynamic and force points 
forces: Aerodynamic forces for each surface 

EQUAtion: Equations of motion coefficient Mi through M 5 ,f^ and^^ 
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Card 3.12.1-Select Frequencies for Printout 

Card 3.12.1 is optional. If omitted, the items requested by card 3.12 will be printed 
for the first reduced frequency only. 


COL. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SELEct 

AlO 

Keyword introducing numbers of frequencies at 
which data will be printed. 

11-15 

IFREQj^ 

1215 

Numbers of frequencies at which data will be 
printed. 

If more than 12 frequencies are desired use two 
SELEct cards. 

16-20 

• 

• 

IFREQ- 
• ^ 

• 



NOTE: The selected array of frequencies applies to all print items. 
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Card 3.12.2-Checkout Printout 

Card 3.12.2 is optional and is used to generate intermediate data for problem 
checkout. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

CHECkout 

AlO 


Keyword requesting the printing of intermediate 
data calculated during the aerodynamics data 
processing . 

11-15 

ICKPRT 

15 

Option indicating the amount of intermediate 
data to be printed. 

=1 Print message when each overlay. is entered 
and indicate field length required for each 
overlay. 

=s3 Print items under 1 above plus many of the 
intermediate calculations 

i.e. 

- for FLXAIC printlAIC] , [AIC] [(|)] r [AICl [(t>" ] 
and the full iFll , lF2] , [fG] . 

- for DLAIC print[CPgg],[CPGgg],[F],[w],[WGl, [ACp] 
and others. 

“ for DLPRES print [Acp] and[Fj^ 2 ] ' 

=5 Print items under 1 and 3 above plus a message 
each time a VARDIM routine is called. 


Cards 3.13 through 3.16 Mode Shape Definition 

The L215 program, modal interpolation (ref. 6) provides all of the mode shapes used 
inside L217 (EOM). The modal data is defined per surface in arrays called {SA} written 
by L215 (ref. 6) onto the magnetic file named ”SATAP.'' 

Cards 3.13 through 3.14.2 are used to define the correspondence between FLEXSTAB 
bodies and the modal interpolation SA arrays. Cards 3.15 and 3.16 serve the same 
purpose for Doublet Lattice bodies/surfaces. It should be noted that order of thin and 
slender bodies is different for the two aerodynamic programs. The correlation of modal 
data and aerodynamic bodies is only concerned with them and slender bodies. 
Interference bodies are not considered, since their effects are accounted for when 
dealing with their respective slender bodies (see sec. 4.0). 
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The thin and slender bodies in FLEXSTAB are ordered in the following manner: 

• Slender bodies on the plane of symmetry. 

• Slender bodies off the plane of symmetry. 

• Thin bodies on the plane of symmetry. 

• Thin bodies off the plane of symmetry. 

Within each catagory of slender and thin bodies, the ordering must match the ordering 
used in the FLEXSTAB aerodynamic geometry program. 

The ordering of the Doublet Lattice data must correspond to the lifting surface 
panel/slender body order established in the Doublet Lattice aerod 3 mamic program 
(L216) which is: 

• Primary vertical lifting surface panel (thin bodies) on the plane of symmetry. 

• Other primary lifting surface panels (thin bodies) off the plane of symmetry. 

• Y-motion slender bodies. 

• Z-motion slender bodies. 

Within each major catagory, the user must assure that the order matches that of the 
Doublet Lattice run which created the aerodynamic data. (Note: When using Doublet 
Lattice aerodynamics, the term "'thin body” refers to a primary lifting surface panel.) 


Cards 3.13 Through 3.14.2-FLEXSTAB Body Modes 

Omit cards 3.13 through 3.14.2 if Doublet Lattice aerodynamics are being used. 

Card 3.13-FLEXSTAB Slender Body Modes 

Repeat card 3.13 for each FLEXSTAB slender body in the order they appear on 
GDTAPE. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-20 

SLENder 

body 

A10,10X 

Keyword introducing a card defining mode shapes 
for a slender body. 

21-25 

ISY 

15 

The numbers of the Modal interpolation Program (ref. 6) surfaces (SA 
arrays) to be used to determine the Y and Z motion, respectively, for 
this slender body. If ISn = 0, modal displacements of zero will be 
used. (See section 4.6.1). 

26-30 

ISZ 

15 

31-35 

ISTORE 

15 

File number (a surface number) of EOMLOD in which 
the aero force matrices will be written. 

ISTORE >0. 

36-40 


5X 

Not Used 

41-50 

XSHIFT 

ElO.O 

Shifting values (defined in the reference axis system) 
which are used to account for the difference between the 
local structural axis system and the local aerodynamic 
axis system, and/or to account for any desired axis 
corrections needed. 

51-60 

YSHIFT 

ElO.O 

61-70 

ZSHIFT 

ElO.O 

j 


♦Normally, ISTORE will be equal to ISY or ISZ. If both ISY and ISZ are defined 
ISTORE should be equal to ISZ. 
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3.14-FLEXSTAB Thin Body Modes 


Repeat cards 3.14 through 3.14.2 for each FLEXSTAB thin body in the order they 
appear on GDTAPE. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-20 

THIN 

body 

A10,10X 

Keyword introducing a card defining mode shapes 
for a thin body. 

21-25 

IS 

15 

Number of the Modal Interpolation Program (ref. 6) surface. (SA 
array ) to be used to determine the motion of this thin body. (See 
section 4.6.1). 

26-30 


5X 

Not Used 

31-35 

ISTORE 

15 

File number (a surface number) of EOMLOD on which the 
aero force matrices will be written.* 

36-40 


5X 

Not Used 

41-50 

XSHIFT 

ElO.O 

Axis shifting values (defined in the reference axis 
system) which account for the difference between the 
local structural and aerodynamic axis system and/or 
any axis correction desired. 

51-60 

YSHIFT 

ElO.O 

61-70 

ZSHIFT 

ElO.O 


*If ISTORE is left blank ISTORE = IS is assumed. If IS = 0 (mode shapes of zero) 
then ISTORE must be non-zero. 
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Cards 3.14.1 and 3.14.2 are optional. They are input in pairs to detine mode shapes 
for subsurfaces (control surfaces) of the previously detined thin body. 

Card 3.14.1-FLEXSTAB Control Surface Modes 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-20 

CONTrol 

surface 

A10,10X 

Keyword introducing a card defining mode shapes 
for a control surface. 

21-25 

IS 

15 

Number of the Modal Interpolation Program (Ref. 6) 
surface corresponding with the motion of the 
control surface. 

26-30 

NCSPAN 

15 

Number of panels on the control surface. 

31-35 

ICSTOR 

15 

j 

File number (a surface number) of EOMLOD in which the 
aero forces will be written. 

If ICSTOR is left blank it is set equal to ISTORE of 
the parent thin body. 

36-70 



Available for conunents. 


Card 3.14.2-Panels on the FLEXSTAB Control Surface 
Card 3.14.2 must follow card 3.14,1. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

IPAN 

1415 

Numbers of the panels of the thin body which are on 

6-10 

IPAN 2 


the control surface. 

IPAN^; i = 1, NCSPAN 

66-70 

IPAN. . 
14 






Repeat the card with 14 numbers per card until all 
panels on the control surface are specified. 


IPANj^ indicates the panel number as counted relative to the first panel on the 
parent thin body. 
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Cards 3,15 Through 3.16-Doublet Lattice Body Modes 


Omit cards 3.15 through 3.16 if FLEXSTAB aerodynamics are being used. 


Card 3.15-Doublet Lattice Thin Body (Primary Lifting Surface Panel) Modes 

Repeat cards 3.15 through 3.15.2 for each Doublet Lattice thin body in the order 
they appear on NGETP. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-20 

THIN 

body 

A10,10X 

Keyword introducing a card defining mode shapes 
for a thin body. 

21-25 

IS 

15 

Number of the Modal Interpolation Program (ref. 6) surface (SA 
array) to be used to determine the motion of this thin body. 

26-30 


5X 

Not Used 

31-35 

ISTORE 

15 

File number (a surface number) of EOMLOD in which 
the aero force matrices will be stored.* 

36-40 


5X 

Not Used 

41-50 

XSHIFT 

ElO.O 

Aerodynamic axis shifting values (defined in reference 

51-60 

YSHIFT 

ElO.O 

axis system) which accounts for the difference between 

61-20 

I 

) 

ZSHIFT 

ElO-O 

the local structural and aerodynamic axis systems and/or 
any axis correction desired. 


* If ISTORE is left blank the program sets ISTORE = IS. 

If IS = 0 (mode shapes of zero) then ISTORE must be non-zero. 


87 




Cards 3.15.1 and 3.15.2 are optional. They are input in pairs; one pair for each 
control surface on the previously defined thin body requiring different mode shape 
definitions. 

Card 3.15.1-Doublet Lattice Control Surface Modes 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

0 

CM 

1 

H 

CONTrol 

surface 

AlO/lOX 

Keyword introducing a card defining mode shapes 
for a control surface. 

21-25 

IS 

15 

Surface number specified to the Modal Interpolation 
Program (Ref. 5) when defining mode shapes for this 
control surface. 

26-30 

NCSBOX 

15 

Number of boxes on the control surface. 

31-35 

ICSTOR 

15 

File number (a surface number) of EOMLOD in which the 
aero forces will be written. 

If ICSTOR is left blank it is set equal to ISTORE of 
the parent thin body. 

36-70 



Available for comments. 


Card 3.15.2-Boxes on a Doublet Lattice Control Surface 
Card 3.15.2 must follow card 3.15.1. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

IBOX^ 

1415 

Numbers of the boxes of the thin body which are on the 
control surface. 

IBOX^ ; i = 1, NCSBOX 

6-10 

IBOX^ 

11-15 

66-70 

IBOXj 

^®°^NCSBOX 


* 

i 

Repeat the card with 14 numbers per card until all 
panels on the control surface have been specified. 


IBOX^ indicates the box number as counted relative to the first 'box on the 
parent thin body. 
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Card 3.16-Doublet Lattice Slender Body Modes 


Repeat card 3.16 for each Doublet Lattice slender body in the order they appear on 
NAETP. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-20 

SLENder 

AlO 

Keyword introducing a card defining mode shapes 


body 

lOX 

for a slender body. 

21-25 

ISY 

15 

. 

Number of the Modal Interpolation Program (ref. 6) surface (SA 
array) to be used to determine the Y or Z motion, respectively. 

26-30 

ISZ 

15 

for this slender body. 

Specify only onel {See section 4.6.1 and 4.7.1) 

31-35 

ISTORE 

15 

File number (a surface number) of EOMLOD in which the 




aero force matrices will be written.* 

36-40 


5X 

Not Used 

41-50 

XSHIFT 

ElO.O 

Aerodynamic axis shifting values (defined in the 

51-60 

YSHIFT 

ElO.O 

reference axis system) which account for the difference 

61-70 

ZSHIFT 

ElO.O 

between the local structural axis and the local aero- 




dynamic axis and/or any axis corrections desired. 

* Normally, ISTORE 

will be 

equal to ISY or ISZ (whichever is defined) . 

Card 3.17-End of Aerodynamic Card Data 

COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

END 

AlO 

Keyword signaling the end of FLEXSTAB or Doublet 




Lattice aerodynamic card input data. 

11-70 



Available for comments 
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Card Set 4.0-Terminate L217 (EOM) Execution 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

0 

1 

$QUIt 

AlO 

Keyword requesting that L217 (EOM) execution be 
terminated . 

11-70 



Available for comments 
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Requirements or Function 

Key Words and/or Variables 

Card Format 

Reference 
Card Set (CS) 


$ROM 

AlO 

1.0 

Job title card (Optional) 

$TITle 

AlO 

1.1 

Job checkout card 
(Optional) 

$NON-stop 

AlO 

1.2 

Begin input of 
Structural Data 

$STRuctural 

AlO 

2.0 

Number of degrees 
of freedom 

SIZE NDOFST NDOFCS 

AlO, 215 

2.1 

Type of analysis 
(Optional 

SYMMetric 
ANT I - s ymme t r i c 

AlO 

2.2 

(Optional) 

PULL 

HALF 

AlO 

2.3 

(Optional) 

CHECkout ICKPRT 

AlO, 15 

2.4 

Modal Interpolation 
! data to use for control 
surface inertial 
; coupling (Optional) 

MODES SATAP 

2A10 

1 

2.5 

Source of Stiffness 
Data 

STIFfness KTAPE KFILE KMAT PRINt 

2A10,2I5,A10 

2.6 j 

Omit card 2.6.1 if KTAPE CARDS (card 2.6) 

Repeat card 2.6.1 for 
each row of K; 
i = 1, NDOFST 

K(i,j) j = 1/ NDOFST 

7E10.0 

2.6.1 

Source of Mass Data 

MASS MTAPE MFILE NMAT PRINt 

2A10,2I5,A10 

2.7 

Omit card 2.7.1 if MTAPE CARDS (card 2.7) 

Repeat card 2.7.1 for 
each row of M 
i = 1, NDOFST 

M(i,j) j = 1/ NDOFS^ 

7E10.0 

2.7.1 


NO 


6.3.4 SUMMARY OF CARD INPUT DATA 




VO 

K> 


Requirements or Function 

Key Words and/or Variable 

Card Format 

Reference 
Card Set (CS) 

Introduce Structural 
Damping (Optional ) 

DAMPing 

AlO 

2.8 

Damping Factors (Must 
follow card 2.8) 

gsd (i) i = l,NDOFST 

7E10.0 

2.8.1 

Omit cards 2.9 and 2.9.1 if NDOFCS = 0 (card 2.1) 

Cards 2.9 and 2.9.1 
are repeated for each 
control surface 
affected by the 
additional control 
surface freedoms. 

SURFace IS NMODEI Tirans Rotl Rot2 I^cal 

A10,2I5,4A10 

2.9 

Repeat card 2.9.1 NNODEI times (card 2.9) 

tlASS^ IMERTY^ INERTX^ 

6E10.0 

2.9.1 

End of Structural Data 

END 

AlO 

2.10 

Introduction of 
Aerodynamic Data 

„ FLEXstab AIC 

DOUBle.t PRESsure 

3 AlO 

I 

3.0 

Aerodynamic constants 

CONStants RHO VEL NTDEF NPDEF 

A10,2E10.0, 

215 

3.1 

Mode Shape Input file 
(Optional) 

SATAp= ISATAP 

2 AlO 

3.2 

FLEXSTAB Aerodynamic 
Data. Read only if 
FLEXSTAB aerodynamic 
option chosen (card 3.0) 
(Optional ) 

GDTApe= GDTAPE IGCASE 

2A10,I5 

3.3 

SDSStp= SDINDX SDDATA ISCASE 

3A10,I5 

3.4 

Doublet Lattice 
Aerodynamic Data. Read 
only if Doublet 
Lattice option chosen 
(card 3. 0) (Optional) 

NGETp= NGETP IGCASF 

2A10,I5 

3.5 j 

NAETp= NAETP lACASE 

2A10,I5 

3.6 { 

1 




Requirements or Function 
(Optional) 

(Optional) 


EOMT ap= 

EOMLod= 


EOMTAP 


EOMLOD 


Quasi-Steady option. 
Only if Doublet Lattice 
AIC Aerodyncimics 
used (card 3.0) 


QUAS i-steady 


Thin Body Control Point 
option. Read only if 
FLEXSTAB Aerodynamics 
used, (card 3.0) 
(Optional) 


CENTr oids 


Key Words and/or Variables 


Reference 

Card Format Card Set (CS) 
2A10 3.7 

2A10 3.8 


AlO 


3.9 


AlO 


3.10 


Omit cards 3.11 and 3.11.1 if Doublet Lattice pressures used 
(card 3.0) 


Gradual Penetration 

GRAD, pen XG NGCP 

AlO.ElO.O, 

15 

3.11 

Gradual Penetration 
Control Points. Repeat 
card 3.11.1 for each 
Control Point i=l, NGCP 

XGCP^ XGLE^ 

2E10.0 

3.11.1 

Print request (Optional) 
Repeat card 3.12 for 
each item desired 

MODES 

PRiNt FORCES 

EQUAtions 

2 AlO 

3.12 

Frequency Selection 
(Optional) 

SELEct IFREQj^ IFREQj . . • 

AlO, 1215 

3.12.1 






Checkout Printout 
(Optional) 


CHECkOUt ICKPRT 


AlO, 15 


3.12.2 





































Requirements or Function 

Key Words and/or Variables 

Card Format 

Reference 
Card Set (CS) 

If Doublet Lattice Aerodynamics are used omit cards 
3.13 through 3.14.2 

Repeat card 3.13 for 
each FLEXSTAB Slender 
Body 

SLENder body ISY ISZ ISTORE XSHIPT YSHIFT 7, SHIFT 

1 

A10,10X,3I5, 

5X,3E10.0 

3.13 

Repeat cards 13.4 
through 3.14.2 for 
each FLEXSTAB Thin Body 

THIN body IS ISTORE XSHIFT YSHIFT ZSHIFT 

A10,10X,I5,5X, 

I5,5X,3E10.0 

3.14 

Control Surface modes 
(Optional) 

CONTrol surface IS NCSPAN ICSTOR 

A10,10X,3I5 

3.14.1 * 

Control Surface panels. 
Must follow card 3.14.1 

IPAN^ i = 1, NCSPAN (card 3.14.1) 

1415 

3.14.2 

I Omit cards 3.15 through 3.16 if FLEXSTAB Aerodynamics are used. ' | 

Repeat cards 3.15 
through 3.15.2 for 
each Doublet Lattice 
Thin Body 

THIN body IS ISTORE XSHIPT YSHIFT ZSHIFT 

A10,10X,I5,5X, 

I5,5X,3E10.0 

3.15 

11 

Control Surface Modes 
(Optional) 

CONTrol surface IS NCSBOX ICSTOR 

A10,10X,3I5 

3.15.1 i' 

Card 3.15.2 must 
follow 3.15.1 

IBOX^ i = 1, NCSBOX (card 3.15.1) 

1415 

3.15.2 

Repeat card 3.16 for 
each Doublet Lattice 
Slender Body 

SLENder body ISY ISZ ISTORE XSHIFT YSHIFT ZSHIFT 

A10,10X,3I5, 

5X,3E10.0 

3.16 

End of Aerodynamic Data 

END 

AlO 

3.17 



$QUIt 

AlO 

4.0 












6.4 MAGNETIC FILE INPUT DATA 


When used as a module in the DYLOFLEX, much of the data required by L217 (EOM) is 
prepared by upstream programs and written on magnetic files (disk or tape) to be read 
and used by L217. The file names and the programs that normally generate them are: 


Default 

Name 

Generating Program 

References 

GSTIFF 

- 

- 

GMASS 

- 

- 

SATAP 

Modal Interpolation L215 

6 

NGETP 

Doublet Lattice L216 

8 

NAETP 

Doublet Lattice L216 

8 

GDTAPE 

FLEXSTAB 

7 

SDINDX 

FLEXSTAB 

7 

SDDATA 

FLEXSTAB 

7 


Note that different file names may be specified via card input data (see cards 2.5 
through 2.7 and 3.2 through 3.6). Also, the files may be generated by different computer 
programs as long as the content and format matches that described on the following 
pages. 

6.4.1 ”GSTIFF’’-GENERALIZED STIFFNESS MATRIX 

"GSTIFF,” the generalized stiffness matrix, may be read by L217 (EOM) from a 
magnetic file (tape or disk) in the READTP/WRTETP format (see ref. 2). L217 (EOM) 
rewinds the file before attempting to read the stiffness matrix [K] . If requested via card 
input data the program will also skip "KFILE” logical files and "KMAT” matrices 
before reading [K] (see fig. 19). 

6.4.2 ^^GMASS”-GENERALIZED MASS MATRIX 

”GMASS,'* the generalized mass matrix, may be read by L217 (EOM) from a magnetic 
file (tape or disk) in the READTP/WRTETP format (see ref. 2). L217 (EOM) rewinds the 
file before attempting to read the mass matrix [M]. If requested via card input data, 
the program will also skip logical "MFILE” files and ”MMAT” matrices before reading 
[M] (see fig. 20). 
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File KFIIE 


File KFTLE+1 


EOF 
[ MATi 3 
[ MAT2 ] 


Matrix KMAT 
Matrix KMAT+1 




[ maTkmat 3 

[K] 



[K] the generalized 
stiffness matrix is 
of size 

NDOFst X NDOFst 


where NDOFst is the 
number of structural 
degrees of freedom. 


Figure 19. — Contents of GSTIFF 



[M] the generalized 
mass matrix is of size 

NDOFst NDOFst 


where 

NDOFST is the number of 
structural degrees of 
freedom. 


Figure 20. — Contents of GMASS 
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6.4.3 "SATAP”-MODAL DATA FILE 


The mode shape data required by L217 (EOM) is generated by the program 
L215 (INTERP) (ref. 6 ) and written on a magnetic file, "SATAP”, (card 2.3) in the 
READTP/WRTETP format described in reference 2. 


L215 (INTERP) writes one logical file of data for each surface for which mode shapes 
are defined. The file number has a one-to- one correspondence with the surface number. 
Thus, when L217 (EOM) needs modal data for surface "IS,” it rewinds "SATAP” and 
skips "IS- 1 ” logical files before reading. 


The contents of one file on "SATAP” are displayed in figure 21. 


Matrix 

Size 

Matrices 
in file i 
of SATAP 

NNOD ES 

X 

NMODES 

C<^^] 

NNODES 

X 

NMODES 

[4>y] 

NNODES 

X 

NMODES 

[♦,] 

NNODES 

X 

NMODES 

t*e ] 

X 

NNODES 

X 

NMODES 

t+e 5 

y 

NNODES 

X 

NMODES 

C*e 1 

Z 

NNODES 

X 

6 

[ qeoro] 

NSA 

X 

1 

(SA) 


End-of-File 


where NMODES = number of modes 

NNODES = number of nodes 

and NSA < 10000 

Note: NSA is dependent on the 

method of interpolation 
(See Reference 5) 


Figure 21. — Contents of a File on SA TAP 


[<^x]> [*^y] [*^z] contain the modal deflections at the nodes in the x, y and z 

Erections of the local structural coordinate system. 

[(^ 0 ^], [<^0 ] and [(^ 0 ^] contain the modal rotations about the x, y and z axes at the 
nodes in the local structural coordinate system. 
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[geom] contains the nodal descriptions. For each node the following six items are 
defined: 

X = coordinate in the reference axis system 
Y = coordinate in the reference axis system 
Z = coordinate in the reference axis system 
</>x = rotation about the x-axis in radians 

<t>y = rotation about the y-axis in radians 
<^2 = rotation about the z-axis in radians 

{SA} contains information to allow interpolation for mode shapes over the nodes. The 
length and contents of {SA} vary with the interpolation type (surface spline, beam 
spline, motion point, motion axis, polynomial). See references for a complete 
description. 

6.4.4 "NGETP”-DOUBLET LATTICE GEOMETRY FILE 

When Doublet Lattice L216 (ref. 8) geometry data is required by L217 (EOM), it is read 
in the READTP/WRTETP format from ''NGETP,” a magnetic file (tape or disk) (see 
ref. 1), L217 (EOM) rewinds the file before attempting to read the data. If the user 
specifies a data case, “IGCASE,” other than number one, the program will also skip 
*TGCASE-1” logical files before reading. 

For each data case, **NGETP” contains 18 matrices followed by an end-of-file (see 
fig. 22). All geometry is defined in the reference axis system. 

The following variables are used to define array sizes on "NGETP.” 

NBOD = number of primary lifting surface panels (thin bodies), plus the number of 
slender bodies 

NBOX = number of aerodynamic boxes on all primary lifting (thin body) and 
interference surfaces 

NTOT = number of aerodynamic elements (boxes plus slender body segments) on all 
surfaces (primary lifting, interference, and slender bodies) 

NSBE = number of slender body segments (counting all slender bodies) 

Three of the matrices on "’NGETP” are described below. The others are self-explanatory. 
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Matrix Size 

Matrix 

Description 

15 X 1 

tGCM} 

Geometry Control Matrix 

8 X NBOD 

[ BODTAB] 

Table of thin and slender bodies 

* X 1 

(PSG} 

Primary surface geometry 

* X 1 

{I SGI 

Interference surface geometry 

♦ X 1 

{SBG} 

Slender, body geometry 

NBOX X 1 

(XCU) 

Coordinates of 1/4 chord 

NBOX X 1 

{YC4} 

points on all thin and 

NBOX X 1 

{ZC4} 

interference surface elements 

NBOX X 1 

{X3C4} 

Coordinates of 3/4 chord 

NBOX X 1 

{Y3C4} 

points on all thin and 

NBOX X 1 

[Z3C4) 

interference surface elements 

NTOT X 1 

{A} 

Areas of all elements 

NTOT X 1 

{Y} 

Dihedral angles of all elements 

NTOT X 1 

1 

{IBSB} 

Table defining correspondence 
between interference body and 
slender body elements 

NS BE X 1 

{SBL} 

Slender body segment lengths 

NSBE X 1 

{SBXl 

Slender body segment midpoint 

NSBE X 1 

{SBR} 

Slender body segment midpoint 
radii 

NSBE X 1 

{SBRS} 

End-of-File 

Slender body segment midpoint 
radius slopes. 


* Difficult to calculate and not required by L217(EOM) 


Figure 22. — Contents of a File on NGETP 


{GCM}-Geometry Control Matrix 


Element 

Type 

Description 

ICASE 

I 

Case number. 

ICOND 

I 

Condition number 

IGEOM 

I 

Geometry data successfully saved if ”IGEOM = 1”. 

REFCHD 

R 

Reference chord length. 

REFSPN 

R 

Reference span length. 

ACAP 

R 

Reference area. 

FMACH 

R 

Mach number. 

NTBB 

I 

Number of thin body boxes. 

NIBB 

I 

Number of interference body boxes. 

NSBE 

I 

Number of slender body segments. 

NTB 

I 

Number of thin bodies (the number of lifting surface 
panels). 

NSB 

I 

Number of slender bodies. 

NPLS 

I 

Number of primary lifting surfaces. 

NPIS 

I 

Number of interference surfaces. 

- 

- 

Not Used. 


where: 

R = Real 

I = Integer 
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[boot AB] -T able of Thin and Slender Bodies 

Note: There is one thin body for each lifting surface panel. The matrix will contain one 
column for each body. 


Element 





Location 

Element 

Type 

Description 

(1,1) 

ITYPE 

I 


1 = slender body and 3 = thin body. 

(2,1) 

IPOS 

I 


0 = off-POS and 1 = on-POS. 

(3,1) 

Xo 

R' 


Origin for local aerodynamic axis 

(4,1) 

Yo 

R 


system of the body (in the referenced 
axis coordinate). 

(5,1) 

Zo 

rJ 



(6,1) 

y 

R 


Dihedral angle of the body. 

(7,1) 

IFIRST 

■1 

1 


Pointers to the first and last boxes 
of the body in the total array of the 

(8,1) 

LAST 

iJ 


structure's boxes. 


where: 


R = Real 

I = Integer 

{IBSB} -Correspondence Between Interference Body Boxes 
and Slender Body Segments 

Each interference body box will be associated with one or two slender body segments: 
two if the physical slender body is represented by two slender bodies, one for y motion 
and another for z motion. Therefore, each element i of the matrix {IBSB} will contain 0, 
1, or 2 numbers. 


If element i represents a thin body box or a slender body segment, then IBSB| = 0. 


If element i represents an interference body box, the first slender body segment with 
which it is associated, JSBl, will appear in the right-most 30 bits of IBSBi. If there is a 
second segment, JSB2, it will appear in the left-most 30 bits. 


IBSBi = 


JSB2 


JSBl 


1 30 


60 
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6.4.5 *‘NAETP”-DOUBLET LATTICE AERODYNAMIC FILE 

When Doublet Lattice L216 (ref. 8) aerodynamic data is required by L217 (EOM), it is 
read from "NAETP,” a magnetic file (tape or disk), in the READTP/WRTETP format 
(see ref. 1). L217 (EOM) rewinds the file before attempting to read the data. If the user 
specifies a data case, "lACASE,” other than number one, the program will skip 
"IACASE-1” logical files on "NAETP” before reading data. 

The number of matrices on “NAETP” for each data case depends upon the option chosen 
in Doublet Lattice. The variable “lAERO” can be extracted from element 3 of the 
aerodynamic control matrix, the first matrix for each case. If: 

lAERO = 1 Matrices [D] and [F] are written onto "NAETP” for each reduced 

frequency 

lAERO = 2 Matrices [ACp] and [Q] are written onto "NAETP” for each 

reduced frequency 

lAERO = 3 Matrices [D], [F], [ACp] and [Q] are all written onto "NAETP” 

The following variables define the size of matrices on "NAETP” displayed in figure 23. 


NFREQM 

NBOX 

NTOT 

NMODE 


= Number of reduced frequencies at which the aerodynamic data is 
defined 

= Number of lifting surface boxes (lifting and interference) 

= Number of elements on all surfaces (lifting, interference, and slender) 
= Number of modes used in Doublet Lattice 
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Matr ix_Sizg 
(r5FFEQ?5+14) X 1 


1 X 2NSBFP 
1 X 2NSBPP 


1 X 2NSBFP 


2NT0T X 1 
2NTOT X 1 


2NT0T X 1 
2NMODE X NMODE 


Ma tr ix 
tACM} 


LfJ 


IPJ; 


IPJ 

I-Matrix 
Trapezoidal 

10HQUASI-END| 

{CP) 1 
{CP} 2 

• 


{CP} 

tQ] 


End-of “File 


Available only if lAERO = 
Available only if lAERO = 
0. nd/or^^are repeated for e 


BSiSElEtion 

Aerodynamic Control Matrix 7se®”next 1 
page) 


[F] matrix by rows 



Two matrices in a variable number of 
records representing the guasi'inverse 

Fecord terminating Quasi“lnverse 

[CP] - pressure coefficient 
Matrix by columns 



Generalized Air Forces 


1 or 3 

2 or 3 

ach reduced frequency; IFREQ =1, NFREQM 


Figure 23. — Contents of a File on NAETP 
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Aerodynamic Control Matrix 


Element 

Location 

Element 

Type 

Description 

1 

ICASE 

I 

Case number. 

2 

ICOND 

I 

Condition number. 

3 

lAERO 

I 

Option indicating the aerodynamic data saved. 
= 1 [D] and [F] saved for each frequency. 

= 2 {CP} and [Q] saved for each frequency. 

= 3 both 1 and 2 above. 

4 

FMACH 

R 

Mach number. 

5 

NFREQM 

I 

Number of reduced frequencies. 

6 

NBOX 

I 

Number of lifting surface boxes. 

7 

NSB 

I 

Number of slender body elements. 

8 

NMODE 

I 

Number of modes. 

9 

REFCHD 

R 

Reference chord. 

10 

REFSPN 

R 

Reference span. 

11 

ACP 

R 

Reference area. 

12 

BR 

R 

Reference length. 

13 

- 

- 

Not used. 

14 

- 

- 

Not used. 

15 

• 

RFREQ(l) 

R 

First reduced frequency. 

# 

• 

RFREQ 

(NFREQM) 

R 

Last reduced frequency. 
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where: 


R = Real 

I = Integer 


6.4.6 ”GDTAPE”-FLEXSTAB GEOMETRY FILE 

Each configuration (i.e., each data case) processed by the GD program is written 
sequentially on "GDTAPE” as one logical file. A file begins with a 50-word header 
record followed by two records for every slender, interference, and thin body in the 
configuration. The file ends with a 50-word trailer record and an end-of-file (EOF) 
mark. A complete description of • **GDTAPE” is provided in appendix A of reference 7, 
volume 2. 

The GD program has the ability to process the geometry definitions for more than one 
airplane configuration; therefore, there may be more than one logical file (i.e., one file 
per aircaft configuration). 

After the GD program terminates, the **GDTAPE” magnetic tape is generated by the job 
control card: 


COPYBF(TAPEl,GDTAPE,n) 

where n = number of logical files generated by the GD program. 

Prior to running L217 (EOM), the data on tape must be recopied to the disk file, named 

"GDTAPE” by the job control card; 

COPYBF(TAPE,GDTAPE,n) 

6.4.7 '•SDINDX/SDDATA”-FLEXSTAB AERODYNAMIC FILE 

SDSSTP is a binary tape generated by the FLEXSTAB SD&SS program (ref. 7) and is 

composed of many logical files. In general, these files consist of; 

• A master catalog of the contents of SDSSTP. This catalog is always the first file on 
the tape. 

• Miscellaneous data required by the TH, SLOADS, and ALOADS programs. These 
data are stored in one file and consist of such information as moments of inertia, 
center of gravity location, number and size of the stability derivative matrices, etc. 

• Slender body and thin body lifting pressure data. These data are stored in one 

logical file. 
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• The [Ki] and [Mx] structural matrices. These matrices (if they exist) are copied 
from SICTP3 onto SDSSTP. 

• Stability derivative matrices. 

• Displacement and camber slope matrices. 

• Elastic axis load matrices (optional). 

• Gust matrices (optional) 

• Any user-specified matrices. An option in the SD&SS program allows the user to 
specify matrices to be stored on SDSSTP, in addition to those that normally are 
saved. 


Each matrix is stored as one logical file. The first file on SDSSTP, is the master catalog. 
It is composed of several data records. There are two t3rpes of data records in this 
catalog-matrix records and terminal records. A matrix record contains information 
pertaining to the data stored in one particular file of SDSSTP. Thus, for each file on 
SDSSTP (not counting the master catalog), there is a corresponding matrix record in the 
master catalog. A terminal record is used to signify the end of the catalog. 

If the SD&SS recycle capability is used (i.e., more than one data case is executed), the 
overall format of SDSSTP changes slightl}". For each successive data case, only those 
matrices that are different from those of the initial case are saved (i.e., no duplication of 
data exists on SDSSTP). This change is also reflected in the first file (the master 
catalog). Additional matrix records corresponding to the files (matrices) added to 
SDSSTP are also included. Terminal records are used to indicate the end of the 
sequence of matrix records for each data case. In other words, the master catalog 
contains: a sequence of matrix records pertaining to the matrices saved from the initial 
case; a terminal record; another sequence of matrix records reflecting the new matrices 
added to SDSSTP from the second data case; another terminal record; and so on, for as 
many cases as executed. A detailed description of SDSSTP is contained in appendix A of 
reference 7, volume 2. 

After execution of the SD&SS program, the SDSSTP tape is generated by the job control 
cards: 


COPYBF(TAPE12,SDSSTP,l) 

COPY(TAPE14,SDSSTP) 

COPYBF(EOF,SDSSTP,l) 


(Master catalog) 

(Matrices, SDSS data, and SPRES data) 
(EOF mark on tape) 
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Before executing L217 (EOM), the data on tape SDSSTP must be copied to disk by the 
job control cards: 

COPYBF(SDSSTP,SDINDX) (Master catalog) 

COPY(SDSSTP,SDDATA) (Matrices required) 

6.5 OUTPUT DATA 

L217 (EOM) will produce a printed output listing and two magnetic files (tape or disk). 

6.5.1 PRINTED OUTPUT 

The program will initialize and terminate the printed output, with special pages 
indicating the program name, version, and date of run. 

All card input data will be printed as it is read. Where appropriate, the card images 
will be augmented by interpretative comments and/or diagnostics. 

The geometry from the aerodynamic program will be printed in tabular form with 
control information based upon card input instructions. 

An explanation of the SIGMA and RG arrays appearing in the geometry tables follows. 

[SIGMA] 

Each row is of [SIGMA] (also known as [S]) indicates where the forces on aerodynamic 
element i are to be stored. The primary purpose of SIGMA is to add the interference 
body forces to the slender body forces. However, as noted later, other factors are 
included in SIGMA for a FLEXSTAB AIC data case. 

SIGMA(i,l) Indicates which force point will receive the SIGMA(i,2) fraction of the 
force of aerodynamic element i. 

SIGMA(i,3) Indicates which force point will receive the SIGMA(i,4) fraction of the 
force of aerodynamic element i. 

For a FLEXSTAB AIC data case, the following ACp factors are included in the second 
and fourth columns of SIGMA: 

• 7t/2 for slender body Y-motion elements. 

• - 7t/2 for slender body Z-motion elements. 

• 1.0 for thin body elements. 

• -1.0 for interference body elements. 
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{RG} 

The gust rotation matrix has elements for all thin body and slender body control points. 
For a thin body element i 

Symmetric analysis: RG = cosy^ 

Antisymmetric analysis; RG = - sin-yj 

where yj is the element's dihedral angle 
For a slender body element i 

Symmetric analysis: Y-motion element RGj = 0. 

Z-motion element RGj = 1. 

Antisymmetric analysis: Y-motion element RGi = -1. 

Z-motion element RGj = 0. 

Note: The gust rotation matrix, which is printed, contains the sign function 6^ shown 
in equation (12). 

All other printed output must be requested with card input options. 

During the processing of structural data, the program will optionally print: 

[K] Generalized stiffness matrix read from tape (see card 2.4). 

[M] Generalized mass matrix read from tape (see card 2.5). 

Cards 3.12 and 3.12.1 may be used to cause the following items to be printed during the 
processing of aerodynamic data: 

• Interpolated mode shapes. 

• Contents of the file ”EOMLOD” (the forces). 

• Contents of the file ’’EOMTAP" (the equations of motion). 

Additional intermediate calculations may be requested using the checkout option, 
card 3.12.2. 
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6.5.2 MAGNETIC FILE OUTPUT 


L217 (EOM) will generate two magnetic files, "EOMTAP” and "EOMLOD.” Both are 
written in the READTPAVRTETP format described in reference 1. 

^'EOMTAP’’ 

The file known as ”EOMTAP” will contain the equations of motion matrices generated 
by the L217 (EOM) program. The matrices are written onto "EOMTAP” by the 
subroutine WRTETP (see ref. 1 for a description of WRTETP formats). The contents of 
"EOMTAP” are displayed in figure 24. 

••EOMLOD” 

The file known as "EOMLOD” will contain aerodynamic geometry and forces for use in 
the DYLOFLEX LOADS program (ref. 7). The data is written onto "EOMLOD” by the 
subroutine WRTETP (see ref. 1 for a description of WRTETP formats). The contents of 
"EOMLOD” are displayed in figure 25. 

Contents of Aerodynamic {HEADER} Array 

The array HEADER will always contain 30 words. 

Word Contents 

1 7HDYLOFLX. 

2 6HL2217VC. 

3 Date of run in format lOHyr/mo/da. 

4 NDOF = number of degrees of freedom (1=^NDOF^100). 

5 NID = 0. 

6 NGCP = number of gust control points (1^NGCP^35). 

7 NFREQM = number of reduced frequencies ^ 1. 

8 q = dynamic pressure = 1/2 

9 V = velocity = true airspeed. 

10-20 Future use = 0. 

21 1 indicates [Mj] is on "EOMTAP.” 

22 0 or 1; a 1 indicates [M 2 ] is on "EOMTAP.” 



size o^_Ma tr ix £iSS£i£&4SQ 

DYLOFLEX Header Matrix describing the 
file 

Generalized stiffness 

Generalized structural damping (Viscous) 
Generalized Mass 

Frequencies. [ Hz ] for which the following 
aerodynamic matrices are defined 

Generalized aerodynamic stiffness 
forces. 

Generalized aerodynamic damping 
forces. 

Distance between gust control point 
and the gust zones 

Excitation aerodynamic generalized 
forces 

[Ms] and f$] are repeated for 
frequencies 2, 3 NFREQM.) 

NDOF = number of generalized coordinates (number of modes) 

NFREQM = number of frequencies at which the aerodynamics are 

defined 

NGCP * number of gust control points (zones) 

1 Contents of HEADER are described on the following page. 

2 [Mj] is present only if requested by Card 2.8. 

3 {fj^} is present only for the first frequency 

4 [ $ ] is complex. Therefore, the actual size written on the file 
ie 2ND0F x NGCP. 

Figure 24. — Contents of EOMTAP 
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size of Matrix 


Matrix 


5®§S£i£tion 


NELE X 2 

[X,Y] 

NELE X 1 

{A1 

NPREQM X 1 

{FREQM} 

NELE X NDOF 

[Fi ] 

NELE X NDOF 

[Fz] 

NELE X NGCP 

[FG] 


End-of-File 

where 



Local coordinates of the surface*s 
elements 

Areas of the surface* s elements 

Frequencies for which the aerodynamic 
matrices below are defined 

Response air forces (coefficients of q) 

Response air forces (coefficients of q) 

Gust air forces (complex) 


[Fi]f [ ^2 ] ] are repeated 

for each frequency 


NELE = 
NDOF = 
NGCP = 


number of aerodynamic elements on surface i. 

number of generalized coordinates 

number of gradual penetration gust zones. 


1 [FG] is complex. Therefore, the actual size written on the file 
is 2N‘ELE X NGCP. 

The logical file shown above is repeated for every surface (see ISTORE 
on Cards 3.12-3.15) • There will be empty files for imbedded surfaces 
not defined. 


Figure 25. — Contents of EOMLOD 



23 1 indicates [M3] is on "EOMTAP.” 

24 1 indicates [M4] is on "EOMTAP.” 

25 1 indicates [M5] is on "EOMTAP.” 

26 0 indicates [Ms] is not on "EOMTAP.” 

27 0 indicates {C2} is not on "ECMTAP.” 

28 0 indicates {C3} is not on "EOMTAP.” 

29 1 indicates {t^} is on "EOMTAP.” 

30 1 indicates [^] is on "EOMTAP.” 

6.6 RESTRICTIONS 


L217 (EOM) problem size maximums are displayed below. 


Variable 

Maximum 

Variable Description 

NMODES 

100 

Number of modes. 

NTOTF 

400 

Total number of aerodynamic elements on all 
bodies (thin, slender, and interference). 

NC3CP 

35 

Number of gust control points (zones). 

NFREQM 

20 

Number of reduced frequencies at which 
aerodynamics are defined. 

Nmax 

100 

' Number of aerodynamic elements for an 
interpolation surface. 


6.7 DIAGNOSTICS 

All errors detected by L217 (EOM) will result in the printing of a diagnostic message. 
The message will be introduced by the line, 

♦♦******FATAL ERROR n IN name 

where "n” indicates the error number (see the list below) and "name” indicates the 
routine printing the message. The introductory line will be followed by .other lines 
describing the error conditions. 
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Error 

Number Description of Error 


1 Premature end-of-file on card input. 

2 Keyword on data card is inappropriate or unrecognizable. 

3 Structural data not processed correctly before processing aerodynamic card 
input data. 

4 Illegal file name. 

5 Illegal file or matrix spacing. 

6 Variable outside allowable range. 

7 FETAD error trying to generate a file’s FIT, buffer, and RA+2 entry. 

8 FETDEI error trying to release a file’s FIT, buffer, and RA-f 2 entry. 

9 READTP error trying to read a matrix from a file (see sec. 6.7.1 for a list of 
READTP error codes.) 

10 WRTETP error trying to write a matrix on a file (see sec. 6.7.2 for a list of 
WRTETP error codes). 

11 STARTR error trying to initialize blank common storage array. 

12 INITIR error trying to generate storage for an array. 

13 DELETE error trying to erase storage of an array. 

15 AINTL or AINTG error trying to interpolate for modes. 

6.7.1 READTP ERROR CODES 
Error Code 

= 0 If no errors are detected during reading. 

= lOOO+I If an FSF (forward space file) error occurred, where I is the number 

of file marks remaining to be skipped when an end-of-information 
was encountered. 

2 If the number of matrices or files to be skipped, before reading 

starts, is less than zero. 
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= 3 

= 3000+1 


= 4 
- 5 


If the dimensioned number of rows in the matrix is less than or 
equal to zero. 

If an FSR (forward space record) error occurred, where I is the 
number of records remaining to be skipped when either an 
end-of«file or end-of-information was encountered. 

Number of rows in the matrix is greater than the dimensioned row 
size in the program. 

If the name check failed. 


= 6 If the number of rows in the matrix (M) times the number of 

columns (N) is greater than the buffer size, or M * N ^ 0. 

= 7 An end-of-file was read. If it occurs while reading the matrix ID, no 

information is stored in the user’s area. If it occurs while reading 
the matrix, the ID information will be stored. Note that the records 
will always be in pairs and an end-of-file should always be 
encountered with the ID record. 

6.7.2 WRTETP ERROR CODES 


Error Code 
= 0 

= lOOO+I 

= 2 
= 3 

= 3000+1 
= 4 


If no errors are detected during writing. 

If an FSF (forward space file) error occurred, where I is the number 
of file marks remaining to be skipped when an end-of-information 
was encountered. 

If the number of matrices or files to be skipped, before writing 
starts, is less than zero. 

If the dimensioned number of rows in the matrix is less than or 
equal to zero. 

If an FSR (forward space record) error occurred, where I is the 
number of records remaining to be skipped when either an 
end-of-file or an end-of-information was encountered. 

If the actual number of rows is greater than the dimensioned 
number of rows in the matrix. 


= 6 


If the number of rows in the matrix (M) times the number of 
columns (N) is greater than the buffer size. 
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7.0 SAMPLE PROBLEM 


The sample problem chosen is the fuselage, wing, and horizontal tail combination shown 
in figure 26. One control surface is placed on the wing. Two rigid body modes (vertical 
trnaslation and pitch) make up the structural degrees of freedom. The control surface 
degree of freedom was added to the basic modes in L215 (INTERP). The FLEXSTAB 
aerodynamic option was chosen for this example. 

Generalized stiffness and mass data is read from file GMASS. Interpolation data is on 
the file MODH7. There are four interpolation arrays. Two rotary inertias are input for 
the control surface due to the rotation of the surface about its swept hinge line. 


Boeing Commercial Airplane Company 
P.O. Box 3707 

Seattle, Washington 98124 
May 1977 
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Figure 26. — Sample Problem Configuration 





Nose UP. 


0 . 

0 . 


WUIT 
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program L217A0 VERSION FEB S, 77 NOW RUNNING. 

• THE PROGRAM IS PART OF THE DVLUFLX SYSTEM 

• DEVELOPED FOR NASA UNDER CONTRACT NAS1>13918. 

‘ DATE OF RUN IS 77/03/02. 

TIME OF RUN IS 18.16.01. 

<*****irt*****t*************************m************* 
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PROGRAM DIRECTIVE CAROS READ BY L217A0 


( $EOM 

lATITLE L2l7(EOMt USAGE DOCUMENT SAMPLE PROBLEM 
(STITLE USING FLEXSTAB AERODYNAMIC DATA FOR H7NC. 

(STITLE MOOES ARE IN THE INERTIAL AXIS SYSTEM. 

I STITLE TRANSLATION IS POSITIVE OOHN> PITCH IS POSITIVE NOSE UP. 
ISSTRUCTURAL DATA 


NOW IN STRUCT TO PREPARE EQUATIONS OF MOTION 
MATRICES (MLf M2t ANO M3J DEPENDENT ONLY UPON 
STRUCTURAL INPUT DATA (MASS* STIFFNESSt ETC.) 


STRUCTURAL CARD INPUT DATA FOLLOWS. 
(SIZE 21 


TOTAL NUMBER OF DEGREES OF FREEDOM =» 3 

( 2) STRUCTURAL AND ( 1) CONTROL SURFACE. 

(SYMMETRIC ANALYSIS 
(HALF AIRPLANE 
(MODES MOOH7 


THE MODE SHAPES WILL BE READ FROM MODH7 


(STIFFNESS GMASS 0 IPRINT 

THE generalized STIFFNESS MATRIX WILL BE READ FROM GMASS 
AFTER SKIPPING 1 MATRICES. 

AND PR INTED AS READ. 


(MASS GMASS 0 OPRINT 


THE GENERALIZED MASS MATRIX WILL BE READ FROM GMASS 
AND PR INTEO AS READ. 


( DAMP ING 


THE damping factors WILL BE READ FROM CARDS 
AND PRINTED AS READ. 


) 


I 

) 

) 


I 


) 


THE MATRIX DAMP 
ROM 1 

( 0. 0. 

I S 

1 BY 

(SURFACE A 

ITRANS 

roty 

CONTROL SURFACE A HAS 1 NODES. 

iranslation freedoms will be used. 

ROTATION DZ/OX FREEDOMS WILL BE USED. 
rotation DZ/OY FREEDOMS WILL BE USED. 
THE COORDINATES ARE IM THE LOCAL AERO 

X Y 

(3.25 4.5 

Z 

0.0 

mass 

0.1 


) 

ROTX LOCAL ) 


XIS SYSTEM. 

INERTIA inertia 

I.O 2.0 I 


( End 


} 
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HATKIX 


LABEL GENERALIZED STIFFNESS 
SIZE 2 ROMS AND 2 COLUMNS 


COLUMN 1 2 

RON 1 0. 0. 


ROW 


2 0 


0 



MATRIX 


LABEL GENERALIZED MASS 

SIZE 2 ROUS AND 2 COLUMNS 


COLUMN 1 2 

BOW I .1763E+03 -.ll00E«-00 

Row 2 -.llOOE^OO .4800E4'01 


STRUCT IS FINISHED. 

MAXIMUM FIELD LENGTH USED BY STRUCT WAS 073343 

RETURN TO L217 DRIVING PROGRAM AND READ ANOTHER PROGRAM DIRECTIVE CARO. 
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PROGRAM OIRECUVE CAROS READ BV L217AO 
( SAERO FLEXSTAB AIC 
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N) 

■U 


NOW IN FLXAIC TO PREPARE EQUATIONS OF MOTION 
BASED UPON FLEXSTAB A IC MATRIX. 


FLEXSTAB CARO INPUT DATA FOLLOWS. 

(CONSTANTS 6.11015E-8 9<9S4.4 1 Z ) 

RHO « .6U02E-07 

VEL » .99&44E+04 

0 - .5*RH0*VEL**2 = .30273E401 


(SATAP MOOH7 > 

MODES WILL BE READ FROM HOOH7 


(GOT APE CDH7 » 

FLEXSTAB GEOMETRY WILL READ FROM GDH7 (CASE 1). 


(SOSSTP SOIH7 SOSH7 I 

FLEXSTAB AEPO DATA WILL BE READ FROM SDIH7 AND SDSH7 (CASE. 1). 


<PRIWT 

MOOES 

I 

< PRINT 

FORCES 

) 

(PRINT 

EQUATIONS 

I 


( SLENDER BODY 2 2 

Y-MOTION WILL NOT BE RETAINED. 

(THIN BODY I .5 1. 

( CONTROL SURFACE 4 1 

( 12 

(THIN BOOY 3 0. 1. 

( END 


0 . 


0 . 


) 


FLEXSTAB CASE IDENTIFICATION - THIS RUN IS A TEST CASE USING THE OOUBLAI TEST CASE H 7WC 

USER IDENTIFICATION * CASE H7WC HAS ONE CONTROL SURFACE 

MACH NUMBER > .8000E<’00 


GO-2 

GO-3 


UNITS ARE IN INCH 



GEOHETRV SUHHARY FOR DOUBLET LATTICE CASE I 
BOX ELEHENT COORDINATES ARE AT QUARTER CHORD. 


syhhetric analysis requested. 


MODAL 


K) 

KJ% 


NO. 

REF. AXIS 
X 

COORDINATES 
Y 1 

DIHEDRAL 
IN DEC. 

AREA 

TYPE 

POS 

SURFACE 
IN OUT 

NO. 

FRACTION NO. 

. FRACTION 

1 

.50 

0.00 

0.00 

0.00 

• 35 

SB 

ON 

0 

2 

1 

•15708£»0l 

0 

0. 

2 

1.50 

0.00 

0.00 

0.00 

• 75 

SO 

ON 

0 

2 

2 

•1570B£*01 

0 

0. 

3 

2.50 

O.QO 

0.00 

0.00 

• 05 

SB 

ON 

0 

2 

3 

.15708E»01 

0 

0. 


3.50 

0.00 

0.00 

0.00 

.93 

SB 

ON 

0 

2 

4 

•15/0UE«01 

0 

0. 

b 

4.50 

0.00 

0.00 

0.00 

.98 

SB 

ON 

0 

2 

5 

•15706E»01 

0 

0. 

6 

5.50 

0.00 

0.00 

0.00 

.95 

SB 

ON 

0 

2 

6 

•1570nE«0l 

0 

0. 

7 

6.50 

0.00 

0.00 

0.00 

• 85 

SB 

On 

0 

2 

7 

•15708C«01 

0 

0. 

6 

7.50 

0.00 

0.00 

0.00 

.75 

SB 

ON 

0 

2 

6 

.15700C*01 

0 

0. 

9 

8.50 

0.00 

0.00 

0.00 

• 60 

SB 

ON 

0 

2 

9 

.l570Bf »0l 

0 0. 

10 

9.50 

0.00 

0.00 

0.00 

.25 

SO 

ON 

0 

2 

10 

•1570BE»01 

0 

0. 

11 

.50 

0.00 

0.00 

0.00 

• 35 

SB 

ON 

2 

2 

11 

-.1570BE»01 

0 

0. 

12 

1.50 

0.00 

0.00 

0.00 

.75 

SB 

ON 

2 

2 

12 

-.l57Gi:£*01 

0 

0. 

13 

2.50 

0.00 

0.00 

0.00 

.05 

SB 

ON 

2 

2 

13 

-.1570UE >01 

0 

0. 

14 

3.50 . 

0.00 

0.00 

0.00 

• 93 

SB 

ON 

2 

2 

14 

-.157U8L*01 

0 

0. 

15 

4.50 

0.00 

0.00 

0.00 

.98 

SB 

ON 

2 

2 

15 

-.l570dE*0l 

0 

0. 

16 

5.50 

0.00 

0.00 

0.00 

.95 

SB 

ON 

2 

2 

16 

-.I5708t»0l 

0 

0. 

17 

6.50 

0.00 

0.00 

0. 00 

.05 

SB 

On 

2 

2 

17 

-.I570dr«0l 

0 

0. 

IB 

7.50 

u.co 

o.uo 

0.00 

.75 

SB 

ON 

2 

2 

la 

-.I5708t«0l 

0 

0. 


8.50 

0.00 

0.00 

0.00 

• 60 

SB 

ON 

2 

2 

19 

-.l570dC*0l 

0 

0. 

20 

9.50 

o.co 

0.00 

0.00 

.25 

SB 

ON 

2 

2 

20 

-.1570QEt01 

0 

0. 

21 

2.50 

.50 

IfOO 

o.no 

1.00 

10 


0 

0 

3 

0. 

13 

-.lOOOOE^Ol 

22, 

3.50 

.50 

1.00 

0.00 

1.00 

1 B 


0 

0 

4 

0. 

14 

-.lOOOOP^Ol 

23 

4.50 

.50 

1.00 

0.00 

1.00 

IB 


0 

0 

5 

0. 

15 

-.10000E*01 

24 

5.50 

• 50 

1.00 

0.00 

1.00 

IB 


0 

0 

6 

0. 

16 

-.lOOOOEfOl 

25 

6.50 

.50 

1.00 

0.00 

1.00 

IB 


0 

0 

7 

0. 

17 

-.lOOOOEtOl 

26 

7.50 

.50 

1.00 

0.00 

1.00 

IB 


0 

0 

B 

0. 

18 

lOOOOE#Ol 

27 

0.5b 

.50 

1.00 

0.00 

1.00 

IB 


0 

u 

9 

0. 

19 

-.lOOOOt^Ol 

20 

9.25 

• So 

1 .00 

0.00 

.50 

IB 


0 

0 

10 

0. 

20 

lOOCOE^Ol 

29 

9.75 

• 50 

1.00 

0.00 

.50 

IQ 


0 

0 

10 

0. 

20 

-. 1000iu*0l 

30 

2.50 

l.OO 

.50 

90.00 

1.00 

IB 


0 

0 

3 

0. 

13 

0. 

31 

3.50 

l.CO 

.50 

90.00 

1.00 

1 B 


0 

0 

4 

0. 

14 

0. 

32 

4.50 

1.00 

.50 

90. 00 

1.00 

IB 


0 

0 

5 

0. 

15 

0. 

33 

5.50 

1.00 

• 50 

90.00 

1.00 

IB 


0 

0 

6 

0. 

16 

0. 

34 

6.50 

l.CO 

.50 

90.00 

1.00 

IB 


0 

0 

7 

0. 

17 

0. 

35 

7.50 

l.CO 

• 50 

90.00 

1.00 

IH 


0 

0 

0 

0. 

IB 

0. 

36 

8.50 

l.CO 

.50 

90.00 

1.00 

lU 


0 

0 

9 

0. 

19 

0. 

37 

9.25 

1.00 

.50 

90.00 

.50 

IB 


0 

0 

10 

0. 

20 

0. 

36 

9.75 

1.00 

.50 

90.00 

.50 

18 


0 

0 

10 

0. 

20 

0. 

39 

‘ 2.50 

1.00 

-.50 

90.00 

1.00 

IB 


0 

0 

3 

0. 

13 

0. 

40 

3.50 

1.00 

-.50 

90. 00 

1.00 

1 e 


0 

0 

4 

0. 

14 

0. 

41 

4.50 

1.00 

-.50 

90.00 

1.00 

IB 


0 

0 

5 

0. 

15 

0. 

42 

5.50 

1.00 

-.50 

90.00 

1.00 

IB 


0 

0 

6 

0. 

16 

0. 

43 

6.50 

l.CO 

-.50 

90.00 

1.00 

1 B 


0 

0 

7 

0. 

17 

0. 

44 

7.50 

1.00 

-.50 

90.00 

1.00 

IB 


0 

0 

8 

0. 

18 

0. 

45. 

e.So 

1.00 

-.50 

90. 00 

1.00 

lU 


0 

0 

9 

0. 

19 

0. 

46 

9.25 

1.06 

-.50 

90. 00 

.50 

1 B 


0 

0 

10 

0. 

20 

0. 

47 

9.75 

1.00 

-.50 

90.00 

.50 

18 


0 

0 

10 

0. 

20 

0. 

48 

2.50 

.50 

-1.00 

160.00 

1.00 

1 8 


0 

0 

3 

0. 

13 

• 10000E40I 

49 

3.50 

.50 

-1 .uo 

180. 00 

1.00 

18 


0 

0 

4 

0. 

14 

• 10000E»01 

50 

4.50 

.50 

-1.00 

160.00 

1.00 

1 U 


0 

0 

5 

0. 

15 

• loocoe »oi 


FAC SF 


OCO O O O O O O OOOOOOOOOOOOOOOOOOOOOM'V^J^Jf^)^4^«lfSlK»(^ir^JrsJ^.>^J^«fu^4^4^»»^« 


SL 

5.50 

.50 

-l.OO 

180.00 

l.OO 

18 

0 

0 

6 

0. 

16 

• 10000E401 

52 

6*50 

• 50 

-l.OO 

180.00 

1.00 

18 

0 

0 

7 

0. 

17 

•lOOOOE 401 

53 

7,50 

.50 

-1.00 

100.00 

1.00 

IB 

0 

0 

8 

0. 

10 

• IOOOOE401 

54 

8.50 

.50 

-l.OO 

100.00 

1.00 

1 B 

0 

0 

9 

0. 

19 

• 100COE401 

55 

9.25 

.50 

-l.OO 

180.00 

.50 

IB 

0 

0 

10 

0. 

20 

• lOCOOF 401 

56 

9.75 

.50 

-l.OO 

1 80. 00 

.50 

I 8 

0 

0 

10 

0. 

20 

•lOOCOCtOl 

57 

4.70 

1.48 

0.00 

0. 00 

.92 

IB OFF 

1 

1 

21 

•10000E»01 

0 

0. 

5fl 

5.62 

1.48 

0-00 

0.00 

.92 

Td OFF 

1 

1 

22 

.lOOOOE^Ol 

0 

0. 

59 

6. 54 

1.48 

0.00 

0.00 

.92 

IB OFF 

1 

1 

23 

.lOOOOE^Ol 

0 

0. 

60 

5.12 

2.48 

0.00 

0,00 

.75 

TB OFF 

1 

1 

24 

.10000E*01 

0 

0. 

CL 

5.87 

2.48 

c.oo 

0.00 

. ?5 

T8 OFF 

1 

1 

25 

• lOOOOE *0 L 

0 

(T. 

62 

6.62 

2.48 

0.00 

0.00 

.75 

TB OFF 

1 

i 

26 

•lOOOOEfOi 

0 

0. 

63 

5.53 

3.48 

0.00 

0.00 

.58 

T» OFF 

1 

1 

27 

.lOOOOE^Ol 

0 

0. 

64 

6.12 

3.40 

0.00 

0. 00 

.58 

TU OFF 

1 

1 

28 

.lOOOOfctOl 

0 

0. 

65 

6.71 

3;46 

0.00 

0.00 

.56 

TB OFF 

1 

I 

29 

• lOOOOE *0L 

0 

0. 

66 

5.94 

4.42 

o.uo 

0.00 

.42 

T0 OFF 

1 

1 

30 

.lOOOOEtOl 

0 

0. 

67 

6.37 

4.47 

0.00 

0.00 

.42 

TB OFF 

1 

1 

3l 

. lOOOOE »0l 

0 

0. 

68 

6.79 

4.47 

0. 00 

0.00 

.42 

TB OFF 

4 

1 

32 

. LOOOOE #0L 

0 

0. 

69 

9.34 

1.48 

0.00 

0. 00 

.44 

TU OFF 

3 

3 

33 

.lOOOOE*Ol 

0 

c. 

70 

9.78 

1.48 

0.00 

0.00 

• 44 

TB OFF 

3 

3 

34 

.lOOOOL «0l 

0 

0. 

71 

9.53 

2.47 

0.00 

0- 00 

• il 

TH OFF 

3 

3 

35 

.lOOOOf »0l 

0 

0. 

72 

9.04 

2.4? 

0.00 

0.00 

.31 

TB OFF 

3 

3 

36 

• lOOOOE «0l 

0 

0. 


rHjfoojrsjrsjisjN^rurvisjoOOOoo 





NO. 

control point G60HETRY 
REF. AXIS COORUINA TE S 
X V Z 

BODY 

type 

POOAL 

INTEKPOLATICN 

SURFACEISI 

KG 

GUST 

ZONE 

1 

.SO 

0. CO 

0.00 

SB 

2 

2 

1.00000 

1 

2 

l.SO 

0.00 

0.00 

so 

2 

2 

1. 00000 

1 

3 

2.50 

O.CO 

0.00 

SB 

2 

2 

1.00000 

1 


3.50 

O.CO 

0.00 

SB 

2 

2 

1.00000 

1 

!> 

6.50 

O.CO 

0.00 

so 

2 

2 

1.00000 

1 

b 

5.50 

O.CO 

0.00 

SB 

2 

2 

l.OOOUO 

1 

7 

6.50 

O.CO 

0.00 

so 

2 

2 

l.OuOOO 

1 

0 

7.50 

O.CO 

0.00 

SB 

2 

2 

1.00000 

1 

9 

8.50 

O.CO 

0-00 

SB 

2 

2 

1.00000 

1 

10 

9.50 

O.CO 

o.uo 

so 

2 

2 

1.00000 

1 

11 

5.02 

1.60 

0.00 

1b 

i 

1 

1.00000 

1 

iZ 

5-96 

1.68 

0.00 

Trt 

1 

1 

1. 00000 

1 

13 

6. 66 

1.68 

0.00 

TU 

1 

1 

I.DUOUO 

1 

U 

5.36 

2.60 

0.00 

TO 

1 

1 

1.00000 

1 

IS 

6.13 

2.68 

0.00 

1U 

1 

1 

1.00000 

1 

16 

6.89 

2.66 

0.00 

TO 

1 

1 

l.OOOUO 

1 

17 

5.76 

3.68 

0.00 

1U 

1 

1 

i.ouono 

1 

IB 

6.32 

3.68 

0.00 

TO 

1 

1 

1.00000 

1 

19 

6.91 

3.68 

0.00 

TH 

1 

L 

1. 00000 

1 

20 

6.09 

6.67 

0.00 

TO 

1 

1 

1.00000 

1 

21 

6.51 

6.67 

0.00 

TB 

1 

1 

1.00000 

1 

22 

6.96 

6.67 

0.00 

TO 

6 

1 

1.00000 

1 

23 

9.69 

1.66 

0.00 

TB 

3 

3 

1.00000 

1 

2S 

9.93 

1.60 

0.00 

TB 

3 

3 

1.00000 

1 

2S 

9.66 

2.67 

0.00 

18 

3 

3 

1.00000 

1 

26 

9.95 

2.67 

0.00 

TU 

3 

3 

1.00000 

1 


/MATRIX 


LABEL MODAL DEFLECTIONS - PARTITION 1 MOOES 
SIZE 26 ROWS AND 3 COLUMNS 


COLUMN 

ROW 

1 

ROW 

2 

ROW 

3 

ROW 

4 

ROW 

5 

ROW 

b 

ROW 

7 

ROW 

8 

ROW 

9 

ROW 

10 

ROW 

1 1 

ROW 

12 

ROW 

13 

ROW 

14 

ROW 

15 

ROW 

16 

ROW 

17 

ROW 

18 

ROW 

19 

ROW 

20 

ROW 

21 

ROW 

22 

ROW 

23 

ROW 

24 

ROW 

25 

ROW 

26 


1 

lOOOE+Ol 
lOOOE + 01 

-.10006+01 
-.10006+01 
lOOOE+01 
-.10006+01 
-.10006+01 
10006+01 
10006+01 
-. 10006+01 
-. lOOOE + 01 
- . 10006+01 
-. 10006+01 
10006 + 01 
10006+01 
-. ICOOE+OI 
-. 10006 +0 I 
-. 100 06+01 
-. 10006+01 
-. 10006+01 
-. 10006 + 01 
-. 10006+01 
-.10006+01 
loooe+oi 

lOOOE+Ol 

-.10006+01 


2 

.A500E+01 

.35006+01 

.25006+01 
.15006+01 
.5 0 0 06 + 00 
-.50006+00 
-.15006+01 
-.25006+01 
-.35006+01 
-.4500E+01 
-.18316-01 
-.9^326+00 

- . 18576+01 
-.3861E+00 
-.11336+01 
-.16866+01 
-.73196+00 
-♦1 3226+01 
-.iSllE+Ol 
- . 10936+01 
-.15156+01 
-.19376+01 
-.44936+01 
-.49346+01 

- .46366 + 01 
-.4953E+01 


0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0. 

0 * 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

-.35406+00 

0 . 

o. 

0 . 

0 . 


1 THROUGH 3 
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I MATRIX 

J' LABEL MODAL SLOPES •> PARTITION 1 MODES 1 THROUGH 3 

SIZE 26 ROWS AND 3 COLUMNS 


COL UMN 

1 

2 

3 

ROW 

1 

0. 

-.IGOOE^-Ol 

0. 

ROW 

2 

0. 

-.10006+01 

0. 

ROW 

3 

0. 

-.10006+01 

0. 

ROW 

4 

0. 

-.lOOOE+Ol 

0. 

ROW 

5 

0. 

-.10006+01 

0. 

ROW 

6 

0. 

-.10006+01 

a. 

ROW 

7 

0. 

-.10006+01 

0. 

ROW 

8 

0. 

-.loooe+oi 

0. 

ROW 

<i 

0. 

-.10006+01 

0. 

ROW 

10 

0. 

-.10006+01 

0. 

ROW 

11 

0. 

-.99926+00 

0. 

ROW 

12 

0. 

- .99986+00 

0. 

ROW 

13 

0. 

-.10066+01 

c. 

ROW 

14 

0. 

-.10006+01 

0. 

ROW 

15 

0. 

-.10006+01 

0. 

ROW 

16 

0. 

-.99916+00 

0. 

ROW 

17 

0. 

- .99666+00 

0. 

ROW 

18 

0. 

- .99676+00 

0. 

ROW 

19 

0. 

-.99676+00 

0. 

ROW 

20 

0* 

-.10016+01 

0. 

Row 

21 

0. 

-.10016+01 

0. 

ROW 

22 

0. 

-.10016+01 

-.98646+00 

ROW 

23 

0. 

-.10006+01 

0. 

ROW 

24 

0. 

-.10006+01 

0. 

ROW 

25 

0. 

-.10006+01 

0. 

ROW 

26 

0. 

-.10006+01 

0. 
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MATRIX 


LABEL FORCE PT . MODAL DEFLECTIONS - PART. 1 MOOES 
SUE 36 ROWS AND 3 COLUMNS 


COLUMN 

X 

2 

3 

ROM 

1 

0. 

0. 

0. 

ROW 

2 

0. 

0. 

0. 

ROW 

3 

0. 

0- 

0. 

ROW 

4 

0. 

0. 

0. 

ROW 

5 

0. 

0. 

0. 

ROW 


0. 

0. 

0. 

ROW 

7 

0. 

0. 

0. 

ROW 

8 

0. 

0. 

0. 

ROW 

9 

0. 

0. 

0. 

ROW 

lO 

0. 

0. 

0. 

ROW 

11 

XOOOEi-01 

•45006+01 

0. 

ROW 

12 

10006 4'Oi 

.3500E+01 

0. 

ROW 

13 

1000E<-01 

.2500E+01 

0. 

ROW 

14 

lOOOE+01 

. 15006 4-01 

0. 

ROW 

15 

ioooe +01 

.5000E+00 

0. 

ROW 

16 

ioooe +01 

-.5000E+00 

0. 

ROW 

17 

1000E4-01 

-.1500E+01 

0. 

ROW 

18 

-• lOOOE^Ol 

-.25006+01 

0. 

ROW 

19 

•.lOOOE+Ol 

-.35006+01 

0. 

ROM 

20 

loooe^oi 

-.45006+01 

0. 

ROW 

21 

10006+01 

.30476+00 

0. 

ROW 

22 

lOoOE+Ol 

- .6209E+00 

0. 

ROW 

23 

lOOOE + Ol 

-.153 86+01 

0. 

ROM 

24 

10006 + 01 

-.12646+00 

0. 

ROW 

25 

-•lOOOE+01 

-.87086+00 

0. 


I Thru 3 
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ROM 

26 

ROW 

27 

ROW 

28 

ROM 

29 

ROW 

30 

ROW 

31 

ROW 

32 

ROW 

33 

Row 

36 

ROW 

35 

POW 

36 


-.lOOOE^Ol 
-.10006+01 
-.10006+01 
-.10006+01 
-.10006+01 
-.10006+01 
-. 10006+01 
-.10006+01 
-.10006+01 
-.10006+01 
-.10006+01 


-.16226+01 

-.32526+00 

-.11166+01 

-.17056+01 

-.9^536+00 

-.13676+01 

-.17896+01 

-.63396+01 

-.67806+01 

-.65256+01 

-.68626+01 


0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

-.20826+00 

0 . 

0 . 

D. 

0. 


TH6 matrices PRIMTEO B6C0W ARE WRITTEM ONTO THE FILE NAMED tOHLOO 



MATRIX —————— — — — — — — w«- 

LAB6L FREQUENCIES AT WHICH THE AERO FORCES ARE DEFINED 
IN CYCLES / SECOND. 

SUE 1 ROUS AND 1 COLUMNS 


COLUMN 1 

0 . 
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MATRIX 


LABEL 

S12E 


ELEMENT LOCAL COORD* (X#YI 
SURFACE 1 

12 ROWS AND 2 COLUMNS 


COLUMN 

ROW 

1 

ROW 

2 

ROW 

3 

ROW 

4 

ROW 

5 

ROW 

6 

ROW 

7 

ROW 

6 

ROW 

9 

ROW 

10 

ROW 

11 

ROW 

12 


1 

• 1202E4-01 
•2121E^01 

• 3040E4'01 
•l617E^0i 
*23706 «^Ol 
•3123E«0l 

• 2032E«^01 

. 26196+01 
*32066+01 

• 2444 ^ 6+01 
•28676+01 
•3289E+01 


2 

•14856+01 
•14656+01 
•14856+01 
•24816+01 
•24816+01 
*24816 + 01 
•3476E+01 

•3476E+01 

*34766+01 

*44676+01 

*44676+01 

*44676+01 



HATRIX 


LABEL 

SIZE 


ELEMENT AREAS 
SURFACE 1 

12 ROUS AND 


COLUMN 

1 

ROW 

1 

• 9167E4-00 

ROW 

2 

•9167E>00 

ROW 

3 

•9167E>00 

ROW 

4 

.7500E4-00 

ROW 

5 

•7500E*00 

ROW 

6 

.7500E4.00 

ROW 

7 

.5833E+00 

ROW 

8 

.5833E+00 

ROW 

9 

•5833E^00 

ROW 

10 

•4167£^00 

ROW 

11 

.4167E+00 

ROW 

12 

.4167E4-00 


I COLUMNS 
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MATRIX 


I 

It 

I 


LABEL 


Fl-RESPONSE 
SURFACE* 1 

AIR FORCES (0 
AT FREQUENCY 1 

COEF.t 

SUE 


12 ROWS 

AND 

i COLUMNS 

COLUMN 


1 


2 

3 

ROW 

1 

0. 


.2865E4-0Z 

.5213E^00 

ROW 

2 

0. 


.9162EV01 

•ASllE^OO 

ROW 

3 

0. 


•ALISE^OI 

.^OSSEt’OO 

ROW 

A 

0. 


•2691E^02 

.63A3E+00 

ROW 

5 

0. 


.7A55E+01 

.5074E4-00 

ROW 

6 

0» 


.3A16E«'0L 

.A981E«00 

ROW 

7 

0. 


.2319E4-02 

.8238E«00 

ROW 

a 

0- 


.6002E+OL 

.7Se9E*00 

Row 

9 

0. 


.2663E«'01 

.991AE4-00 

ROW 

10 

0. 


.1T36E4'02 

.L169E»01 

ROW 

11 

0. 


.3835E«-01 

.1S3AE4-01 

ROW 

12 

0. 


.lESAE't-Ol 

.5895E>01 
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MAHUX 


LABEL F2-RESP0NSE AIR FORCES (Q DOT COEF.) 
SURFACE. 1 AT FREQUENCY I 


SIZE 


12 ROMS AND 

COLUMN 

1 

ROU 

1 

.28796-02 

ROW 

z 

.91878-03 

ROW 

3 

.4065E-03 

RON 

4 

.27038-02 

ROW 

5 

.74866-03 

ROW 

6 

.36306-03 

ROW 

7 

.23326-02 

ROW 

8 

.6029E-03 

ROW 

9 

.28756-03 

ROW 

10 

. 1744E-02 

ROW 

11 

.3652E-03 

ROW 

12 

.16616-03 


3 COLUMNS 


2 

3 

.1873E-02 

.18796-04 

.2796E-02 

•16266-04 

.22^9E-02 

.14556-04 

.2^456-02 

.22876-04 

.2156E-02 

.18296-04 

.1678E-02 

.17966-04 

.2649E-02 

.29706-04 

.1569E-02 

.27366-04 

•1136E-02 

.35746-04 

.23406-02 

.42156-04 

.90416-03 

.55306-04 

.60866-03 

.21256-03 
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HATRIX 


LABEL FGUST- GUST AIR FORCES ICONPLEXI 
SURFACEf I AT FREQUENCV 1 

SIZE 12 ROMS AND 1 COLUMNS 

I COMPLEX) 


COLUHN 



1 

AON 

L 

•2879E-02 

0. 

ROW 

2 

•9187E-03 

0. 

ROW 

3 

•406SE-03 

0. 

ROW 

4 

.2703E-02 

0. 

ROW 

5 

.7486E-03 

0. 

ROW 

6 

.3630E-03 

0. 

ROW 

7 

.2332E-02 

0. 

ROW 

t 

.6029E-03 

0. 

ROW 

9 

.2875E-03 

0. 

ROW 

lO 

.17446-02 

0. 

ROW 

11 

.38526-03 

0. 

ROW 

12 

. 1661E-03 

0. 



MATRl 

LABEL 

SUE 

X — — 

CLEMENT LOCAL COORO. 
SURFACE 2 

10 ROWS ANO 2 

(XtVI 

COLUMNS 

COLUMN 

1 

2 


(tow 

I 

.5000E«^00 

0. 


ROW 

2 

•ISOOE^Ol 

0. 


ROW 

3 

.2500E+01 

0« 


ROW 


.3500E*01 

0. 


ROW 

5 

.ASOOE^Ol 

0. 


ROW 

6 

.55OOE»0l 

0. 


ROW 

7 

.6500E+01 

0. 


ROW 

8 

.TSOOe^Ol 

0. 


ROW 

9 

.8$00E»0l 

0. 


ROW 

10 

,9500e»0l 

0. 
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MATRIX 


LABEL 

SUE 


ELEMENT AREAS 
SURFACE 2 


10 ROWS AND 


COLUMN 


1 

ROW 

1 

•3SOOE^OO 

ROW 

2 

• 7SC0E «>00 

RON 

3 

• 85006 1-00 

ROW 

A 

•9250E^00 

ROW 

S 

•9750£^00 

ROW 

6 

•9500E»00 

ROW 

7 

•eSOOE^OO 

ROW 

8 

.TSOOEtOO 

ROW 

9 

•bOOOE^OO 

ROW 

10 

•2500E^C0 


i COLUMNS 



MATRIX 


LABEL 


fi-rcsponsr a ir forces IQ 

SURFACE.t 2 AT FREQUENCY 

COEF.I 

1 

sreE 


20 ROWS 

AND 3 COLUMNS 

COLUMN 

1 

2 

3 

ROW 

1 

0. 

0. 

0 . 

ROM 

2 

0. 

0. 

0* 

RON 

3 

0« 

0. 

0. 

RON 

A 

0. 

0« 

0. 

RON 

5 

0. 

0. 

0. 

ROW 

6 

0. 

0* 

0. 

ROW 

7 

0. 

0. 

0. 

ROW 

8 

0. 

0. 

0. 

ROW 

9 

0. 

0* 

0. 

ROW 

10 

0. 

0. 

0. 

Ron 

11 

0. 

.5733E+01 

0. 

ROW 

12 

0. 

.5529E+01 

0. 

ROW 

13 

0. 

.7lllE»0l 

.2481E«’00 

ROW 

14 

0. 

.6840e«01 

.2782£»00 

ROW 

15 

0. 


.5107E»00 

Ron 

16 

0. 

,9193E«-0l 

.6590E«00 

ROW 

17 

0. 

• 9291E4-01 

.6l9ie«-00 

ROW 

18 

0. 

.S811E«00 

.9368E.«00 

Aon 

19 

0. 

-.2807E«’Ol 

.3383E«00 

ROW 

20 

0, 

-.2«05E«-0l 

.2621E«00 
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HATHIX 


LABEL F2-RESP0NSE AlA FORCES fQ OOT COEF.I 
SUKFACEt Z AT FREQUENCY I 

SIZE 20 RONS AND 3 COLUMNS 


COLUMN 

1 

2 

3 

RON 

i 

0. 

0. 

0. 

RON 

2 

0. 

0. 

0. 

RON 

3 

0. 

0. 

0. 

RON 

4 

0. 

0. 

0. 

ROW 

5 

0, 

0. 

0. 

ROW 

6 

0. 

0. 

0. 

RON 

7 

0. 

0. 

0. 

ROW 

8 

0. 

0. 

0. 

ROW 

9 

0. 

0. 

0. 

ROW 

10 

0. 

0. 

0. 

ROW 

11 

.5760E-03 

-.2041E-02 

0. 

ROW 

12 

•5554E-03 

-.1119E-02 

0. 

ROW 

13 

.71426-03 

•1331E-02 

.8944E-05 

ROW 

14 

. 8879E-03 

•1S53E-02 

.1003E-04 

ROW 

15 

.1272E-02 

.3145E-02 

.1841E-04 

ROW 

16 

.9224E-03 

.3475E-02 

.23766-04 

Row 

17 

.4298E-03 

.2662E-02 

•22326-04 

ROW 

id 

.580 IE- 04 

.1 5356-02 

.15756-04 

ROW 

19 

-.2812E-03 

.43836-03 

.12206-04 

ROW 

20 

-.2405E-03 

.24356-03 

.94516-05 
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MATRIX 


LABEL PGUST- cusr AIR FORCES ICONPLEX) 
SURFACEt 2 AT FREQUENCY 1 

SIZE 20 ROWS AND I COLUMNS 

(COMPLEX) 


COLUMN 

AOW 

1 

I 

0. 

1 

0. 

ROW 

2 

0. 

0. 

ROW 

3 

0. 

0. 

ROW 

4 

0. 

0. 

ROW 

5 

0. 

0. 

ROW 

6 

0. 

0. 

ROW 

7 

0. 

0. 

ROW 

8 

0. 

0. 

ROW 

9 

0. 

0. 

ROW 

10 

0. 

0. 

ROW 

11 

.5760E-03 

0. 

ROW 

12 

.55546-03 

0. 

ROW 

13 

.7142E-03 

0. 

ROW 

19 

.88796-03 

0. 

ROW 

15 

.12726-02 

0. 

ROW 

16 

.9224E-03 

0. 

ROW 

17 

.42986-03 

0. 

ROW 

18 

.58016-04 

0. 

ROW 

19 

-.29126-03 

0. 

ROW 

20 

-.24056-03 

0. 
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MATRIX 


label element local coord. (X«VI 

SURFACE 3 

SIZE 4 ROUS AND 2 COLUMNS 


COLUMN 


1 

2 

RON 

I 

•3393E^OO 

• 1476E4^0| 

ROM 

2 

•7798E*00 

• L47CiE4>01 

ROW 

3 

.525OE^00 

•2467E«01 

ROW 

4 

.04l7e^OO 

• 2467E4>01 



HATR IX 


LABEL 

ELEMENT AREAS 
SURFACE 3 

SUE 

A ROWS ANO 

COLUMN 

1 

RON 

1 .A3T5E^OO 

ROW 

2 .A375C+00 

ROW 

3 .3125E«00 

ROW 

4 .312SE«'00 


1 COLUMNS 
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HATRIX 

LABEL 

sue 


COLUMN 


1 

2 

3 

ROW 

1 

0. 

•5A28£«00 

•43646^00 

ROW 

2 

0. 

•lOOOE^OO 

•1093E^00 

ROM 

3 

0. 

•1402E^0l 

•4224E^OO 

ROM 

4 

0. 

• 2i87E>00 

•8212E-0L 


fl-R£SPONSE AIR FORCES (0 
SUKfACEt 3 AT FREQUENCY 


COEFO 


4 ROMS AND 


3 COLUMNS 



MATRIX 


LABEL F2-RESP0NSE AIR FORCES IQ DOT COEF.) 
surface, 3 AT FREOLENCV 1 

SIZE A ROMS AND 3 COLUMNS 


COLUMN 


1 

2 

3 

ROM 

I 

.6043E-0A 

•4C59E*02 

.15736-04 

ROM 

2 

.10a9E~0A 

•1154E-02 

•39466-05 

RON 

3 

. IA09E-03 

•35046-02 

.15236-04 

ROW 

4 

.2219E-OA 

.62146-03 

.29616-05 
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MATRIX 


LABEL FGUST- GUST AIR FORCES (COMPLEX! 
SURFACE# 3 AT FREQUENCY I 


SIZE 4 ROWS AND 1 COLUMNS 

(COMPLEX) 


COLUMN 



L 

ROW 

1 

.6043E-04 

0. 

ROW 

2 

•1089E-04 

0. 

ROW 

3 

• 1409E-03 

c. 

RON 

4 

.2219E-04 

0. 


THE MATRICES PRINTED BELOW ARE WRITTEN ONTO THE FILE NAMED EOMTAP 
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MATRIX 


LABEL Ml 

SIZE 3 ROWS AMO 3 COLUMNS 


COLUMN I 2 3 

ROW 1 0. 0. 0. 

ROW 2 0 • 0« 0« 

ROW 3 0. 0. 0. 
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HATRIX 


LABEL 

SIZE 

M? 

3 Rows 

AND 

3 COLUMNS 

COLUMN 


I 

2 

ROk/ 

1 

0. 

0. 

0, 

ROW 

2 

0. 

0. 

0 

ROW 

3 

0. 

0. 

0 
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MATAfX 


LA6EL M3 

SUE 3 ROWS AND 3 COLUMNS 


COLUMN 
ROW 1 

ROW 2 

ROW 3 


1 

.3526E+03 

-. 22006+00 

.16CAE-01 


2 . 

-.2200E+00 

.9600E+01 

.10146+01 


3 

1644E-01 

10146+01 

10306+01 
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MATRIX 

LABEL FREQM 

SIZE 1 ROWS AND 1 COLUMNS 


1 


COLUMN 


0 


MATRIX 


LABEL 

SIZE 

M4 

i ROWS ANO 

3 COLUMNS 


COLUMN 


i 

2 

3 


1 

0. 

.322584^03 

.3383E402 

AOW 

2 

0. 

.6A07E+02 

,5046E»02 

ROW 

3 

0. 

.6887(4^00 

.2455E401 
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MATRIX 


LASEi. N5 

SIZE 3 ROWS AMO 3 COLUMNS 


tOtUMN 


1 

2 

3 

ROW 

1 

•3239E-01 

.7551E-01 

.1220E-02 

ROW 

2 

.6482E-02 

.IA17E*00 

.1619E-02 

ROW 

3 

.^918E~04 

.2S35E-03 

.8851E-04 


MATRIX 


LABeL DISTANCES FROM GUST REF. PT. TO TONES. 
SUE 1 ROWS AND 1 COLUMNS 


COLUMN 1 

0 . 
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MATRIX 


LABEL PHIT 

SIZE 3 ROWS AND 1 COLUMNS 

I COMPLEX ) 


COLUMN 


1 

ROW 

1 

-.3239E-01 

0. 

ROW 

2 

-,^482E-02 

0. 

ROW 

3 

-.6918E-04 

0. 


maximum *F lELO ^length USED BY FLXAIC WAS 072551 
RETURN TO L217 TO READ ANOTHER DIRECTIVE CARD. 
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P«OGRAM DIRECTIVE CAROS READ BY I.217AO 
( iOUlT 
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